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Preface 


This manual describes the routines in release 2.0 of the Image Processing Library, which have been 
developed to facilitate image processing on AMT DAP Series machines. The routines are mainly 
written in FORTRAN-PLUS, and users are assumed to be familiar with this language. 


The principles of FORTRAN-PLUS are outlined in: 

DAP Series: Introduction to FORTRAN-PLUS (man001) 

Detailed information on writing FORTRAN-PLUS programs is supplied in the AMT reference 
publication: — 


DAP Series: FORTRAN-PLUS Language . : (man002) 


Developing and running programs on the DAP is described in the following publications: 
DAP Series: Program Development under UNIX (man003) 


DAP Series: Program Development under VAX/VMX (man004) 


This is the third edition of the Image Processing Library. 
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Chapter 1 


Introduction 


1.1 Overview 


The Image Processing Library has been designed to facilitate image processing on the AMT DAP. 
The DAP (standing for ‘Distributed Array of Processors’) is a massively parallel computer which 
attaches to a host as a peripheral processor. 


The processors in the DAP are arranged in a square matrix; for example 32 x 32 in the case of the 
DAP 500 range and 64 x 64 in the DAP 600. The number of processors on one side of the square 
is called the edge-size of the DAP. ES is used as an abbreviation for edge-size in this manual. 


This manual describes the routines included in the AMT Image Processing Library. The routines 
are divided into four chapters: Image Conversion Routines, Image Processing Primitives, Low 
Level Image Processing Routines and Image Analysis Routines, a summary of their contents is. 
listed below. Howto access the library is described in section 1.2. This is followed by a list of the 
routines in alphabetical order in chapter 2. 


Chapter 3. Image Conversion Routines 


There are three standard ways of storing images in DAP store: raster, sheet and crinkled mappings. 


e By raster images we mean that the pixels of the image have been assigned a row-major 
ordering. Raster images are mapped into the DAP store by imposing long-vector (column- 
major) ordering on each DAP matrix and then column-major ordering on the set of matrices. 
Each pixel in the image is stored in the corresponding position in DAP store, the first pixel 
in the first position etc.. 


e A sheet mapped image is divided into tiles — each tile fills one DAP matrix. A set of tiles is 
stored in column-major ordering in DAP memory. 


e A crinkle mapped image is conceptually divided into ESx EStiles. The size of each tile is: 


height of the image : width of the image 
ES ES 
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Each tile is then ‘crinkled’ up and stored in column-major order under the corresponding 
processor element. Thus each DAP matrix holds one pixel value from each tile. 


Image conversion routines convert between these three ways of mapping the image in DAP store. 
The image must be of size 2“ x2" where 5 < N,M < 10, and may be any (positive) number of 
bits deep, restricted only by storage considerations. Most of the routines in the following chapters 
assume that the image is held in DAP store using a crinkled mapping. Conversions between two’s 
complement and unsigned formats are also included, since calculations in the library routines are 
performed using two’s complement arithmetic. 


CRINK. RASTER — converts an image from crinkled to raster ordering in DAP store. 
CRINK_SHEET — converts an image from crinkled to sheet ordering in DAP store. 
RASTER. CRINK — converts an image from raster to crinkled ordering in DAP store. 
RASTER. SHEET — converts an image from raster to sheet ordering in DAP store. 

SHEET. CRINK — converts an image from sheet to crinkled ordering in DAP store. 
SHEET. RASTER ~— converts an image from sheet to raster ordering in DAP store. 

TWO.UNSIG — converts an image from from two’s complement to unsigned integer 


format maintaining the ordering of the grey-scale values. 


UNSIG.TWO- ~ converts an image from from unsigned to two’s complement integer 
format maintaining the ordering of the grey-scale values. 


Chapter 4 Image Processing Primitives 


Image processing primitives are routines for adding, subtracting, and shifting entire images. Rou- 
tines specified as 8-bits per pixel routines may also be used for 16-bits per pixel images by replacing 
the 8 in the routine name by 16 (and 16 if present by 32). For example: 


MULT.8.TO_16 — multiplies two 8-bit two’s complement images to form a 16-bit 
image. 

MULT-16. TO 32 ~ multiplies two 16-bit two’s complement images to form a 32-bit 
image. 


The routines ADD.8, SCAL.ADD-8 and SUB.8 return their results modulo 256. The routines 
ADD.16, SCAL.ADD-_16 and SUB.16 return their results modulo 65536. 


Arithmetic operations 


ADD.8 - adds together two 8-bit two’s complement images. The result is 
returned modulo 256. 
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Chapter 4 Image Processing Primitives (continued) 


SCAL.ADD.8 7 


SUB_8 = 


MULT.8.TO_16 _ 
SCAL.MULT_8.TO-16 - 


SCAL.DIV_8 - 


Shifting operations 


adds an 8-bit two’s complement scalar to an 8-bit two’s comple- 
ment image. The result is returned modulo 256. 


subtracts one 8-bit two’s complement image from another. The 
result is returned modulo 256. 


multiplies two 8-bit two’s complement images to form a 16-bit 
image. 


multiplies an 8-bit two’s complement image by an 8-bit two’s com- 
plement scalar to form a 16-bit image. 


divides an 8-bit two’s complement image by an 8-bit two’s com- 
plement scalar to form an 8-bit image. 


The shifting routines are for crinkled mapped images only. Shifts north are documented here but 
shifts to other directions can be made by replacing NORTH in the routine name by SOUTH, EAST 


or WEST as required. 


SHIFT.IMAGE.NORTH.P - 


SHIFT_.IMAGE.NORTH.C - 


SHIFT.ROW.NORTH.P - 


SHIFT.ROW.NORTH.C - 


SHIFT.COL.NORTH.P  - 


SHIFT.COL.NORTH.C - 


SHIFT.SHEET.NORTH_P - 


SHIFT.SHEET.NORTH.C —- 
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shifts an entire image north with planar boundary conditions. Ze- 
roes are introduced at the boundaries. 


shifts an entire image north with cyclic (wrap-around) boundary 
conditions. 


shifts a row of matrices of an image north with planar boundary 
conditions. Zeroes are introduced at the boundaries. 


shifts a row of matrices of an image north with cyclic (wrap- 
around) boundary conditions. 


shifts a column of matrices of an image north with planar bound- 
ary conditions. Zeroes are introduced at the boundaries. 


shifts a column of matrices of an image north with cyclic (wrap- 
around) boundary conditions. 


shifts a matrix of an image north with planar boundary conditions. 
Zeroes are introduced at the boundaries. 


shifts a matrix of an image north with cyclic (wrap-around) 
boundary conditions. 
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Chapter 5 Low Level Image Processing Routines 


Low-level image processing routines perform tasks such as edge detection routines, convolutions 
and Fast Fourier Transforms (FFTs). Routines specified as 8-bits per pixel may be used for 16-bits 
per pixel images by replacing the 8 in the name by 16, except for the COMP.FFT_2D routines 
which perform FFTs on images with 8, 16 and 24 bits per pixel. 


ABS. THRESH. 8 _ 
AVERAGE. 8 _ 


BOXLIN. BOX_8 7 


CO6_LFFTLESS = 


CO6.FFT.LV _ 


COMP_FFT_2D_8.TO.8 


COMP_FFT._2D.8.TO_16 


COMP._FFT_2D_8.TO.24 


COMP. FFT. 2D_16.TO-_16 - 


COMP-_FFT_2D_16_TO-.24 - 


thresholds an 8-bit two’s complement image. 


convolves a square averaging mask with an 8-bit two’s complement 
image. 


applies box-in-a-box prescreening to an 8-bit per pixel two’s com- 
plement image. The result is compared with a user supplied 
threshold and a binary result returned. 


calculates the two dimensional discrete Fourier transform of 
ES x ES complex points. 


performs a one dimensional Fast Fourier Transform of ES? com- 
plex points. 


calculates the two dimensional fast Fourier transform of an 8- 
bit two’s complement complex image producing an 8-bit two’s 


. complement complex image. The input data may be scaled; the 


output transform is scaled automatically. 


calculates the two dimensional fast Fourier transform of an 8- 
bit two’s complement complex image producing a 16-bit two’s 
complement complex image. The input data may be scaled; the 
output transform is scaled automatically. 


calculates the two dimensional fast Fourier transform of an 8 
bit two’s complement complex image producing a 24-bit two’s 
complement complex image. The input data may be scaled; the 
output transform is scaled automatically. 


calculates the two dimensional fast Fourier transform of a 16- 
bit two’s complement complex image producing a 16-bit two’s 
complement complex image. The input data may be scaled; the 
output transform is scaled automatically. 


calculates the two dimensional fast Fourier transform of a 16- 
bit two’s complement complex image producing a 24-bit two’s 
complement complex image. The input data may be scaled; the 
output transform is scaled automatically. 
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Chapter 5 Low Level Image Processing Routines (continued) 


COMP_FFT_2D_24. TO_24 — 


COMP_FFT_2D.REAL.3 - 


CONVOLVE._8 = 


DIFF.OF.GAUSS.8 - 


F01.G.MM _ 


FO01. M_INV _ 


F04.QR_GIVENS.SOLVE - 
FILL_IN-1 = 


HISTOGRAM. C_8 _ 
HISTOGRAM.S._8 - 


KIRSCH_8 - 


LAPLACE.8 ~ 


LINE. DET.8 - 


NORMALIZE. 8 a 
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calculates the two dimensional fast Fourier transform of a 24- 
bit two’s complement complex image producing a 24-bit two’s 
complement complex image. The input data may be scaled; the 
output transform is scaled automatically. 


calculates the two dimensional fast Fourier transform of a 24-bit 
floating point complex image producing a 24-bit floating point 
complex image. The input data may be scaled; the output trans- 
form is scaled automatically. 


calculates the convolution of an 8-bit per pixel two’s complement 
image with a square mask. The result may be scaled to avoid 
overflow. 


convolves one of a set of nine difference-of-Gaussian masks with 
an 8-bit two’s complement image. 


performs a general matrix multiply of two matrices A and B where 
A is a PxQ matrix and B is a QxR matrix with P,Q and R in 
the range 1 to ES. 


calculates, in place, the inverse of a given Nx N matrix with N in 
the range 1 to ES. 


solves the linear system Az = ) for z where A is an nxn matrix 
with 2<n<ES+1. The routine may be used to simultaneously 
solve for up to E'S different right hand side vectors b. 


performs a region growing process on a logical image starting from 
input seed points. The region growing 1s constrained by the input 
image — it can only grow where the image is .TRUE.. 


calculates the histogram for an 8-bit two’s complement crinkle 
mapped image or tiles thereof. 


calculates the histogram for an 8-bit two’s complement sheet 
mapped image or tiles thereof. 


applies a user-selected Kirsch compass gradient mask to an 8-bit 
per pixel two’s complement image. The result may be scaled to 
avoid overflow. 


applies a user-selected Laplacian mask to an 8-bit per pixel two’s 
complement image. The result may be scaled to avoid overflow. 


applies a user-selected line detection mask to an 8-bit per pixel 
two’s complement image. The result may be scaled to avoid over- 
flow. 


expands the range of grey-scale values used within an 8-bit two’s 
complement image to the full range. 
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Chapter 5 Low Level Image Processing Routines (continued) 


PERC. THRESH. 8 = 


PREWITT.8 = 


PSEUDO. MEDIAN. 8 - 


PURE.MEDIAN.8  - = 


ROBERTS.8 - 


SOBEL-. 8 _ 


ZERO. X_8 ; — 


Chapter 6 Image 


performs thresholding on an 8-bit two’s complement image. 


applies a user-selected Prewitt compass gradient mask to an 8-bit 
per pixel two’s complement image. The result may be scaled to 
avoid overflow. 


calculates the pseudo-median at each pixel in an 8-bit two’s com- 
plement image for a rectangular neighbourhood centred on the 
image. 


calculates the true median at each pixel in an 8-bit two’s com- 
plement image for a rectangular neighbourhood centred on the 
image. 


applies a user-selected Roberts edge detection mask to an 8-bit 
per pixel two’s complement image. The result may be scaled to 
avoid overflow. 


sums the convolution of an 8-bit per pixel two’s complement image 
with the 3x3 Sobel vertical and horizontal edge detection masks. 
The result may be scaled to avoid overflow. 


finds the zero crossings in an 8-bit two’s complement image. 


Analysis Routines 


Image analysis routines allow the interesting features on an image to found and identified auto- 
matically. SEGMENT_16 and FEATURES. 16 are available for analysing 16-bit images. 


SEGMENT_8 - 


LABEL. 16 _ 


FEATURES_8 = 


CLASSIF = 


extracts the most ‘interesting’ blob or related set of blobs from an 
8-bit two’s complement image by successively applying difference- 
of-Gaussian convolution operators to the image. 


labels the distinct ‘blobs’ in a binary image, returning a 16-bit 
two’s complement image with all the pixels in each blob having 
one value, different for each blob. 


calculates invariant moments and other features for use in classifi- 
cation of ‘blobs’ of each blob in an 8-bit two’s complement labelled 
image (see LABEL- 16). 


calculates the class expectation index (how likely a ‘blob’ is to be 
in a given class) for each set of input blob features. 
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1.2 Access to the library 


Routines in the image processing library are linked in at the consolidation stage of the compiling 
process. For details of compiling and linking see DAP Series: Program Development Under UNIX, 
or DAP Series: Program Development Under VAX/VMS. 


1.2.1 Using the library under UNIX 


Routines are linked in by specifying iplib with the -1 (for library) flag in either dapa or dapf. For 
example: 


dapf -o myfile.dd myfile.df -1 iplib 


dapf will compile the FORTRAN-PLUS program myfile.df, linking in any routines from the 
image processing library which are required and produce a DOF file myfile.dd. dapf and dapa 
will automatically find the version of the library to match the size of DAP that the program is 
being produced for. 


1.2.2 Using the library under VAX/VMS 


You can link routines from the image processing library into a program by using the /LIBRARY 
qualifier to the DLINK command. Two versions of the image processing library are supplied, 
IPLIB5 for DAP 500, and IPLIB6 for DAP 600; when linking routines from the image processing 
library into your program you need to specify the appropriate version of IPLIB. 


For example, to compile and link a FORTRAN-PLUS program in the file PICTURE. DFP for a 
DAP 600 you can use the following commands: 


$ DFORTRAN/DAPSIZE=64 PICTURE 
$ DLINK/DAPSIZE=64 PICTURE, SYS$LIBRARY : IPLIB6/LIBRARY 


To compile and link the FORTRAN-PLUS program in file IMAGES.DFP for a DAP 500 the 
commands would be: 


$ DFORTRAN/DAPSIZE=32 IMAGES 
$ DLINK/DAPSIZE=32 IMAGES ,SYS$LIBRARY : IPLIB5/LIBRARY 


You can define the logical name DAPn-LIBRARY by using the command: 
$ DEFINE DAPn.LIBRARY SYS$LIBRARY:IPLIBn 


where n is 5 (for DAP 500) or 6 (for DAP 600). This will cause DLINK to search IPLIBn 
automatically for unsatisfied external references. If you are going to use IPLIBn frequently, you 
could insert the DEFINE command above into your LOGIN.COM file. If there are several DAP 
users on the system, linked to a DAP 600 say, the system manager could include the command: 
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$ DEFINE/SYSTEM DAP6_LIBRARY SYS$LIBRARY: IPLIB6.DLB 


into the site system start-up command file which would give all users automatic access to the 
library. : 


Similarly, the command: 
$ DEFINE/SYSTEM DAPS5_LIBRARY SYS$LIBRARY:IPLIB5.DLB 


would do the same thing for a DAP 500 system. 


On asystem that has both DAP 500 and DAP 600, then both DAP5_LIBRARY and DAP6_LIBRARY 
can be defined, and users would pick up the version of IPLIB to match the /DAPSIZE that had 
been specified. 
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Chapter 2 


Alphabetical listing of routines 


Name 
ABS. THRESH-_8 ° 
ADD_8 


AVERAGE._8 


BOX_IN. BOX_8 


CO6_FFT_ESS 


C06. FFT.LV 


CLASSIF 


COMP_FFT.2D.8.TO.8 


COMP. FFT_2D_8_TO-16 


Image Processing Library 


Function 
thresholds an 8-bit two’s complement image. 


adds together two 8-bit two’s complement images. The result 
is returned modulo 256. 


convolves a square averaging mask with an 8-bit two’s comple- 
ment image. 


applies box-in-a-box prescreening to an 8-bit per pixel two’s 
complement image. The result is compared with a user sup- 
plied threshold and a binary result returned. 


calculates the two dimensional discrete Fourier transform of 
ES x ES complex points. 


performs a one dimensional Fast Fourier Transform of ES? 
complex points. 


calculates the class expectation index (how likely a ‘blob’ is to 
be in a given class) for each set of input blob features. 


calculates the two dimensional fast Fourier transform of an 8- 
bit two’s complement complex image producing an 8-bit two’s 
complement complex image. The input data may be scaled; 
the output transform is scaled automatically. 


calculates the two dimensional fast Fourier transform of an 8- 
bit two’s complement complex image producing a 16-bit two’s 
complement complex image. The input data may be scaled; 
the output transform is scaled automatically. 
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74 
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Alphabetical listing of routines 





Name Function Page 


COMP. FFT_2D_8.TO_24 -—- calculates the two dimensional fast Fourier transform of an 8- 80 
bit two’s complement complex image producing a 24-bit two’s 
complement complex image. The input data may be scaled; 
the output transform is scaled automatically. 


COMP_FFT_2D_.16_.TO_16 — calculates the two dimensional fast Fourier transform of a 16- 83 
bit two’s complement complex image producing a 16-bit two’s 
complement complex image. The input data may be scaled; 
the output transform is scaled automatically. 


COMP. FFT_2D_16.TO-_24 - calculates the two dimensional fast Fourier transform of a16- 86 
bit two’s complement complex image producing a 24-bit two’s 
complement complex image. The input data may be scaled; 
the output transform is scaled automatically. 


COMP_FFT.2D_24_TO.24 - calculates the two dimensional fast Fourier transform of a24- 89 
bit two’s complement complex image producing a 24-bit two’s 
complement complex image. The input data may be scaled; 
the output transform is scaled automatically. 


COMP_FFT.2D_-REAL.3 - calculates the two dimensional fast Fourier transform ofa 24 92 
bit floating point complex image producing a 24-bit floating 
point complex image. The input data may be scaled; the out- 
put transform is scaled automatically. 





CONVOLVE.8 - calculates the convolution of an 8-bit per pixel two’s comple- 94 
ment image with a square mask. The result may be scaled to 
avoid overflow. 


CRINK. RASTER ~ converts an image from crinkled to raster ordering in DAP 16 
store. 

CRINK. SHEET - converts an image from crinkled to sheet ordering in DAP 18 
store. 

DIFF.OF.GAUSS_8 - convolves one of a set of nine difference-of-Gaussian masks with 96 


an 8-bit two’s complement image. 


F01.G.MM - performs a general matrix multiply of two matrices A and B98 
where A is a PxQ matrix and B is a QxR matrix with P,Q 
and R in the range 1 to ES. 


F01.M.INV - calculates, in place, the inverse of a given Nx N matrix with N 100 
in the range 1 to ES. 


F04.QR.GIVENS.SOLVE -— solves the linear system Az = b for z where Ais annxnmatrix 102 
with 2<n<ES +1. The routine may be used to simultaneously 
solve for up to E'S different right hand side vectors 0. 
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Alphabetical listing of routines 





Name Function Page 


FEATURES_8 — calculates invariant moments and other features for use inclas- 138 
sification of ‘blobs’ of each blob in an 8-bit two’s complement 
labelled image (see LABEL- 16). 


FILL-IN.1 — performs a region growing process on a logical image starting 104 
from input seed points. The region growing is constrained by 
the input image — it can only grow where the image is .TRUE.. 


HISTOGRAM. C_8 — calculates the histogram for an 8-bit two’s complement crinkle 106 
mapped image or tiles thereof. 


HISTOGRAM.S._8 — calculates the histogram for a 8-bit two’s complement sheet 108 
mapped image or tiles thereof. 


KIRSCH_8 — applies a user-selected Kirsch compass gradient mask to an 110 
8-bit per pixel two’s complement image. The result may be 
scaled to avoid overflow. 


LABEL. 16 — labels the distinct ‘blobs’ in a binary image, returning a 16-bit 136 
two’s complement image with all the pixels in each blob having 
one value, different for each blob. 





LAPLACE_8 ~— applies a user-selected Laplacian mask to an 8-bit per pixel 112 
two’s complement image. The result may be scaled to avoid 
overflow. 

LINE. DET.8 — applies a user-selected line detection mask to an 8-bit per pixel 114 
two’s complement image. The result may be scaled to avoid 
overflow. 

MULT_8. TO. 16 — multiplies two 8-bit two’s complement images to form a 16-bit 40 
image. 

NORMALIZE. 8 — expands the range of grey-scale values used within an 8-bit 116 
two’s complement image to the full range. 

PERC_THRESH_8 — performs thresholding on an 8-bit two’s complement image. 118 

PREWITT-_8 — applies a user-selected Prewitt compass gradient mask to an 120 


8-bit per pixel two’s complement image. The result may be 
scaled to avoid overflow. 


PSEUDO. MEDIAN_8 — calculates the pseudo-median at each pixel in an 8-bit two’s 122 
complement image for a rectangular neighbourhood centred 
on the image. 


PURE. MEDIAN.8 — calculates the true median at each pixel in an 8-bit two’s com- 124 
plement image for a rectangular neighbourhood centred on the 
image. 
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Name 
RASTER. CRINK _ 


RASTER. SHEET _ 
ROBERTS-8 = 
SCAL_ADD-_8 _ 
SCAL.DIV_8 - 
SCAL.MULT.8. TO. 16 - 


SEGMENT-_8 = 


SHEET. CRINK - 


SHEET. RASTER ~ 
SHIFT.COL.NORTH.C  - 


SHIFT.COL.NORTH.P  - 
SHIFT_IMAGE.NORTHLC - 


SHIFT.IMAGE.NORTH_P - 


SHIFT.ROW_.NORTH.C - 
SHIFT.ROW.NORTH.P - 
SHIFT.SHEET.NORTH-.C -— 


SHIFT.SHEET.NORTH-.P - 


SOBEL 8 ~ 
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Alphabetical listing of routines 


Function 


converts an image from raster to crinkled ordering in DAP 
store. 


converts an image from raster to sheet ordering in DAP store. 


applies a user-selected Roberts edge detection mask to an 8-bit 
per pixel two’s complement image. The result may be scaled 
to avoid overflow. 


adds an 8-bit two’s complement scalar to an 8-bit two’s com- 
plement image. The result is returned modulo 256. 


divides an 8-bit two’s complement image by an 8-bit two’s 
complement scalar to form an 8-bit image. 


multiplies an 8-bit two’s complement image by an 8-bit two’s 
complement scalar to form a 16-bit image. 


extracts the most ‘interesting’ blob or related set of blobs 
from an 8-bit two’s complement image by successively applying 
difference-of-Gaussian convolution operators to the image. 


converts an image from sheet to crinkled ordering in DAP 
store. 


converts an image from sheet to raster ordering in DAP store. 


shifts a column of matrices of an image north with cyclic (wrap- 
around) boundary conditions. 


shifts a column of matrices of an image north with planar 
boundary conditions. Zeroes are introduced at the boundaries. 


shifts an entire image north with cyclic (wrap-around) bound- 
ary conditions. 


shifts an entire image north with planar boundary conditions. 
Zeroes are introduced at the boundaries. 


shifts a row of matrices of an image north with cyclic (wrap- 
around) boundary conditions. 


shifts a row of matrices of an image north with planar bound- 
ary conditions. Zeroes are introduced at the boundaries. 


shifts a matrix of an image north with cyclic (wrap-around) 
boundary conditions. 


shifts a matrix of an image north with planar boundary con- 
ditions. Zeroes are introduced at the boundaries. 


sums the convolution of an 8-bit per pixel two’s complement 
image with the 3x3 Sobel vertical and horizontal edge detec- 
tion masks. The result may be scaled to avoid overflow. 


Page 


20 


22 


126 


36 


44 


42 


134 


24 


26 


56 


04 


48 


46 


O2 


50 


60 


58 


128 
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Alphabetical listing of routines 





Name Function . Page 


SUB. 8 ~ subtracts one 8-bit two’s complement image from another. The 38 
result is returned modulo 256. : 


TWO.UNSIG - converts an image from from two’s complement to unsigned in- 28 
teger format maintaining the ordering of the grey-scale values. 


UNSIG. TWO — converts an image from from unsigned to two’s complement in- 30 
teger format maintaining the ordering of the grey-scale values. 


ZERO_X_8 - finds the zero crossings in an 8-bit two’s complement image. 130 
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Chapter 3 


Image Conversion Routines 


There are three standard ways of storing images in DAP store: raster, sheet and crinkled mappings. 


e By raster images we mean that the pixels of the image have been assigned a row-major 
ordering. Raster images are mapped into the DAP store by imposing long-vector (column- 
major) ordering on each DAP matrix and then column-major ordering on the set of matrices. 
Each pixel in the image is stored in the corresponding position in DAP store, the first pixel 
in the first position etc.. . 


e A sheet mapped image is divided into tiles - each E'S x ES pixel tile is stored in one DAP 
matrix. A set of tiles is stored in column-major ordering in DAP memory. 


e A crinkle mapped image is conceptually divided into ES? tiles of size (height /E.S )x (width /ES). 
Each tile is then ‘crinkled’ up and stored in column-major order under the corresponding 
processing element. Thus each DAP matrix holds one pixel value from each tile. 


Image conversion routines convert between these three ways of mapping the image in DAP store. 
The image must be of size 2@ x2" where 5 < N,M < 10, and may be any (positive) number of 
bits deep, restricted only by storage considerations. Most of the routines in the following chapters 
assume that the image is held in DAP store using a crinkled mapping. Conversions between two’s 
complement and unsigned formats are also included, since calculations in the library routines are 
performed using two’s complement arithmetic. 
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CRINK. RASTER Image Conversion Routines 


3.1 CRINK_RASTER 


16 


Purpose 


This routine converts an image from crinkled to raster ordering in DAP store. 


Specification 
CRINK. RASTER (IMAGE, BITS , NS_SIZE , WE_SIZE) 
LOGICAL IMAGE(,, BITS , NS_SIZE , WE_SIZE) 


INTEGER *4 BITS, NS_SIZE, WE_SIZE 


Description 
This routine converts images (in place) from crinkled to raster ordering in DAP store where: 


crinkled mapping divides up an image into E'S? tiles of size (height /ES ) x (width /ES ) 
each of which is then ‘crinkled up’ and stored under one processing element. Thus each 
DAP matrix in which the image is stored holds one pixel value from each tile. The number 
of bits per pixel is specified in BITS. Here DAP matrix refers to BITS consecutive bit 
planes 


raster images are mapped into DAP store by imposing long-vector (column-major) or- 
dering on each DAP matrix of pixel values and then column-major ordering on the set 
of matrices 


References 


None. 


Arguments 
IMAGE 
IMAGE is the image to be remapped into DAP store. 


BITS 
On entry BITS is the number of bits per pixel in IMAGE. Unchanged on exit. 


NS.SIZE 


On entry NS_SIZE contains the height of the image in pixels divided by ES. NS_SIZE 
must be one of 1, 2, 4, 8, 16, 32. Unchanged on exit. 


WE.SIZE 


On entry WE-SIZE contains the width of the image divided by ES. WE.SIZE must be 
one of 1, 2, 4, 8, 16, 32. Unchanged on exit. 


Errors 


No explicit error checking is done. 


Auxiliary Routines 
LOG2, X05LOG2 
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Image Conversion Routines CRINK. RASTER 


8 Accuracy 
Not applicable. 


9 Further Comments 
This routine uses Parallel Data Transforms. 


10 Keywords 
Parallel Data Transforms 
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CRINK. SHEET Image Conversion Routines 


3.2 CRINK_SHEET 
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Purpose 


This routine converts an image from crinkled to sheet ordering in DAP store. 
Specification 

CRINK. SHEET (IMAGE, BITS , NS_SIZE , WE_ SIZE) 

LOGICAL IMAGE(,, BITS , NS_SIZE , WE_SIZE) 


INTEGER *4 BITS, NS_SIZE, WE_SIZE 


Description 
This routine converts images (in place) from crinkled to sheet ordering in DAP store where: 


crinkled mapping divides up an image into ES? tiles of size (height /HS ) x (width /ES ) 
each tile is then ‘crinkled up’ and stored under one processing element. Thus each DAP 
matrix in which the image is stored’holds one pixel value from each tile. The number 
of bits per pixel is specified in BITS. Here DAP matrix refers to BITS consecutive bit 
planes 


sheet mapping divides the image up into tiles of ES x ES pixels, each of which occupies 
one DAP matrix. The set of DAP matrices are stored in column-major order 


References 


None. 


Arguments 
IMAGE 
IMAGE is the image to be remapped into DAP store. 
BITS 
On entry BITS is the number of bits per pixel in IMAGE. Unchanged on exit. 


NS_SIZE 


On entry NS_SIZE contains the height of the image in pixels divided by ES. NS_SIZE 
must be one of 1, 2, 4, 8, 16, 32. Unchanged on exit. 


WE-SIZE 


On entry WE-SIZE contains the width of the image divided by ES. WE_SIZE must be 
one of 1, 2, 4, 8, 16, 32. Unchanged on exit. 


Errors 


No explicit error checking is done. 


Auxiliary Routines 
LOG2, X05LOG2 
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8 Accuracy 
Not applicable. 


9 Further Comments 


This routine uses the Parallel Data Transforms. 


10 Keywords 


Parallel Data Transforms 
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RASTER. CRINK Image Conversion Routines 


3.3 RASTER_CRINK 
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Purpose 


This routine converts an image from raster to crinkled ordering in DAP store. 
Specification 

RASTER. CRINK (IMAGE, BITS , NS_SIZE , WE_SIZE) 

LOGICAL IMAGE (,, BITS ,NS_SIZE 'WE.SIZE) 


INTEGER *4 BITS, NS.SIZE, WE.SIZE 


Description 
This routine converts images (in place) from raster to crinkled ordering in DAP store where: 


raster images are mapped into DAP store by imposing long-vector (column-major) or- 
dering on each DAP matrix of pixel values and then column-major ordering on the set 
of matrices. Here DAP matrix refers to BITS consecutive bit planes 


crinkled mapping divides up an image into E'S? tiles of size (height /ES )x (width /ES ) 
each tile is then ‘crinkled up’ and stored under one processing element. Thus each DAP 
matrix in which the image is stored holds one pixel value from each tile. The number of 
bits per pixel is specified in BITS 


References 


None. 


Arguments 
IMAGE 
IMAGE is the image to be remapped into DAP store. 
BITS 
On entry BITS is the number of bits per pixel in IMAGE. Unchanged on exit. 


NS_SIZE 


On entry NS-SIZE contains the height of the image in pixels divided by ES. NS_SIZE 
must be one of 1, 2, 4, 8, 16, 32. Unchanged on exit. 


WE-SIZE 


On entry WE-SIZE contains the width of the image divided by ES. WE_SIZE must be 
one of 1, 2, 4, 8, 16, 32. Unchanged on exit. 


Errors 


No explicit error checking is done. 


Auxiliary Routines 
LOG2, X05LOG2 
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8 Accuracy 
Not applicable. 


9 Further Comments 


This routine uses Parallel Data Transforms. 


10 Keywords 
Parallel Data Transforms 
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RASTER. SHEET Image Conversion Routines 


3.4 RASTER_SHEET 
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Purpose 


This routine converts an image from raster to sheet ordering in DAP store. 
Specification 

RASTER. SHEET (IMAGE, BITS , NS_SIZE , WE_SIZE) 

LOGICAL IMAGE(,, BITS , NS_SIZE , WE_SIZE) 


INTEGER *4 BITS, NS_SIZE, WE.SIZE 


Description 
This routine converts images (in place) from raster to sheet ordering in DAP store where: 


raster images are mapped into DAP store by imposing long-vector (column-major) or- 
dering on each DAP matrix of pixel values and then column-major ordering on the set 
of matrices. Here DAP matrix refers to BITS consecutive bit planes 


sheet mapping divides the image up into tiles of ES x ES pixels, each tile occupies one 
DAP matrix. The set of DAP matrices are stored in column-major order 


References 


None. 


Arguments 
IMAGE 
IMAGE is the image to be remapped into DAP store. 


BITS 
On entry BITS is the number of bits per pixel in IMAGE. Unchanged on exit. 


NS.SIZE 


On entry NS_SIZE contains the height of the image in pixels divided by ES. NS_SIZE 
must be one of 1, 2, 4, 8, 16, 32. Unchanged on exit. 


WE. SIZE 


On entry WE_SIZE contains the width of the image in pixels divided by HS. WE.SIZE 
must be one of 1, 2, 4, 8, 16, 32. Unchanged on exit. 


Errors 


No explicit error checking is done. 


Auxiliary Routines 
LOG2, X05LOG2 
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8 Accuracy 
Not applicable. 


9 Further Comments 


This routine uses Parallel Data Transforms. 


10 Keywords 
Parallel Data Transforms 
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SHEET. CRINK | Image Conversion Routines 


3.56 SHEET.CRINK 
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Purpose 


This routine converts an image from sheet to crinkled ordering in DAP store. 
Specification 

SHEET. CRINK (IMAGE, BITS , NS_SIZE , WE. SIZE) 

LOGICAL IMAGE(,, BITS, NS_SIZE , WE_SIZE) 


INTEGER *4 BITS, NS.SIZE, WE_SIZE 


Description 
This routine converts images (in place) from sheet to crinkled ordering in DAP store where: 


sheet mapping divides the image up into tiles of ES x ES pixels, each of which occupies 
one DAP matrix. The set of DAP matrices are stored in column-major order. Here DAP 
_matrix refers to BITS consecutive bit planes 


crinkled mapping divides up an image into ES? tiles of size (height /ES ) x (width /ES ) 
each of which is then ‘crinkled up’ and stored under one processing element. Thus matrix 
in which the image is stored holds one pixel value from each tile. The number of bits per 
pixel is specified in BITS 


References 


None. 


Arguments 
IMAGE 
IMAGE is the image to be remapped into DAP store. 
BITS 
On entry BITS is the number of bits per pixel in IMAGE. Unchanged on exit. 
NS_SIZE 


On entry NS_SIZE contains the height of the image in pixels divided by ES. NS.SIZE 
must be one of 1, 2, 4, 8, 16, 32. Unchanged on exit. 


WE_SIZE 


On entry WE-.SIZE contains the width of the image divided by ES. WE.SIZE must be 
one of 1, 2, 4, 8, 16, 32. Unchanged on exit. 


Errors 


No explicit error checking is done. 


Auxiliary Routines 
LOG2, X05LOG2 
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8 Accuracy 
Not applicable. 


9 Further Comments 


This routine uses the Parallel Data Transforms. 


10 Keywords 
Parallel Data Transforms 
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SHEET. RASTER Image Conversion Routines 


3.6 SHEET.RASTER 
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Purpose 


This routine converts an image from sheet to raster ordering in DAP store. 


Specification 


SHEET. RASTER (IMAGE, BITS , NS_SIZE , WE.SIZE) 
LOGICAL IMAGE(,, BITS, NS_SIZE , WE.SIZE) 


INTEGER *4 BITS, NS_SIZE, WE_SIZE 


Description 
This routine converts images (in place) from sheet to raster ordering in DAP store where: 


sheet mapping divides the image up into tiles of ES x ES pixels, each tile occupies one 
DAP matrix. The set of DAP matrices are stored in column-major order. Here DAP 
matrix refers to BITS consecutive bit planes 


raster images are mapped into DAP store by imposing long-vector (column-major) or- 
dering on each DAP matrix of pixel values and then column-major ordering on the set 
of matrices 


References 


None. 


Arguments 
IMAGE 
IMAGE is the image to be remapped into DAP store. 
BITS 
On entry BITS is the number of bits per pixel in IMAGE. Unchanged on exit. 
NS.SIZE 


On entry NS_SIZE contains the height of the image in pixels divided by ES. NS_SIZE 
must be one of 1, 2, 4, 8, 16, 32. Unchanged on exit. 


WE. SIZE 


On entry WE-.SIZE contains the width of the image divided by FS. WE.SIZE must be 
one of 1, 2, 4, 8, 16, 32. Unchanged on exit. 


Errors 


No explicit error checking is done. 


Auxiliary Routines 
LOG2, X05LOG2 


Accuracy 
Not applicable. 
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9 Further Comments 


This routine uses the Parallel Data Transforms. 


10 Keywords 


Data Transformation 
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TWO. UNSIG Image Conversion Routines 


3.7 TWO_UNSIG 
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Purpose 


This routine converts an image from from two’s complement to unsigned integer format 
maintaining the ordering of the grey-scale values. 


Specification 
TWO-_UNSIG (IMAGE, BITS , NS_SIZE , WE. SIZE) 
LOGICAL IMAGE(,, BITS,NS_SIZE,WE_SIZE) 


INTEGER *4 BITS, NS_SIZE, WE_SIZE 


Description 


This routine converts images from two’s complement to unsigned format by flipping the most 
significant bit of the pixel value. 


References 


None. 


Arguments 
IMAGE 
IMAGE is the image to be format converted. 
BITS 
On entry BITS is the number of bits per pixel in IMAGE. Unchanged on exit. 
NS_SIZE | 


On entry NS_SIZE contains the height of the image in pixels divided by ES. Unchanged 
on exit. 


WE. SIZE . | 
On entry WE.SIZE contains the width of the image divided by ES. Unchanged on exit. 


Errors 


No explicit error checking is done. 


Auxiliary Routines 


None. 


Accuracy 
Not applicable. 


Further Comments 


None. 
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10 Keywords 


Image Processing 


Image Processing Library 
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UNSIG. TWO Image Conversion Routines 


3.8 UNSIG.TWO 
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Purpose 


This routine converts an image from from unsigned integer format to two’s complement 
integer format maintaining the ordering of the grey-scale values. 


Specification 
UNSIG. TWO (IMAGE, BITS,NS_SIZE, WE-_ SIZE) 
LOGICAL IMAGE (, , BITS,NS_SIZE,WE-SIZE) 


INTEGER *4 BITS, NS_SIZE, WE. SIZE 


Description 


This routine converts images from unsigned to two’s complement format by flipping the most 
significant bit of the pixel value. 


References 


None. 


Arguments 
IMAGE 
IMAGE is the image to be format converted. 
BITS 
On entry BITS is the number of bits per pixel in IMAGE. Unchanged on exit. 
NS_SIZE 


On entry NS_SIZE contains the height of the image in pixels divided by ES. Unchanged 
on exit. 


WE. SIZE 
On entry WE_SIZE contains the width of the image divided by ES. Unchanged on exit. 


Errors 


No explicit error checking is done. 


Auxiliary Routines 


None. 


Accuracy 
Not applicable. 


Further Comments 


None. 
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10 Keywords 


Image Processing 
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Chapter 4 


Image Processing Primitives 


These routines provide image processing primitives such as image shifting and addition of two 
images. Routines specified as 8-bits per pixel routines may also be used for 16-bits per pixel 
images by replacing the 8 in the routine name by 16 (and 16 if present by 32). For example: 


MULT_8.TO_16 - multiplies two 8-bit two’s complement images to form a 16-bit 
image. 

MULT_16. TO. 32 - multiplies two 16-bit two’s complement images to form a 32-bit 
image. 





The routines ADD.8, SCAL.ADD-_8 and SUB.8 return their results modulo 256. The routines 
ADD_16, SCAL.ADD_16 and SUB_16 return their results modulo 65536. 


The shifting routines are for crinkled mapped images only. Shifts north are documented here but 
shifts to other directions can be made by replacing NORTH in the routine name by SOUTH, EAST 
or WEST as required. 
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ADD_8 Image Processing Primitives 


Arithmetic Operations 


4.1 ADD_8 
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Purpose 


This routine adds together two 8-bit two’s complement images. The result is returned modulo 
256. 


Specification 
ADD-.8 (IMAGE1,IMAGE2,NS_SIZE,WE_SIZE,IMAGE.SUM) 


INTEGER*1 IMAGE1 (, ,NS.SIZE,WE- SIZE), 
& IMAGE2 (, ,NS_SIZE,WE-SIZE), 
& IMAGE.SUM (,,NS_SIZE,WE_SIZE) 


INTEGER*4 NS_SIZE,WE- SIZE 


Description 
IMAGE1 is added modulo 256 to IMAGE2 pixel-by-pixel to form IMAGE._.SUM 


References 


None. 


Arguments 
IMAGE1 


On entry IMAGEI is one of the 8-bit images to be added together. Each dimension of 
the image size must be > ES and a multiple of ES. IMAGE1 may be overwritten by the 
sum if desired. 


IMAGE2 


On entry IMAGE2 is the other image to be added. The dimensions of IMAGE2 must be 
identical to those of IMAGE1. IMAGE2 may be overwritten by the sum if desired. 


NS_SIZE 


On entry NS_SIZE contains the height in pixels of the images divided by ES. Unchanged 
on exit. 


WE.SIZE 


On entry WE-SIZE contains the width in pixels of the image divided by ES. Unchanged 
on exit. 


IMAGE_SUM 
On exit IMAGE.SUM is the sum of IMAGE1 and IMAGE2 modulo 256. 


Errors 


No explicit error checking is done. 
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7 Auxiliary Routines 


None. 


8 Accuracy 
Not applicable. 


9 Further Comments 


None. 


10 Keywords 


Image Processing 
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4.2 SCAL_.ADD_8 
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Purpose 


This routine adds an 8-bit two’s complement scalar to an 8-bit two’s complement image. 
The result is returned modulo 256. 


Specification 
SCAL. ADD-_8 (IMAGE,SCAL,NS_SIZE,WE-.SIZE IMAGE_SUM) 


INTEGER*1 IMAGE (,, NS_SIZE,WE-SIZE), 
& -  §CAL, 
& IMAGE.SUM (,, NS_SIZE,WE-SIZE) 


INTEGER*4 NS_SIZE,WE- SIZE 


Description 
SCAL is added modulo 256 to each pixel of IMAGE to form IMAGE_SUM. 


References 


None. 


Arguments 


IMAGE 


On entry IMAGE is the 8-bit image to be added to a scalar. Each dimension of the 
image size must be >ESand a multiple of ES. IMAGE may be overwritten by the sum 
if desired. 


NS.SIZE 


On entry NS_SIZE contains the height in pixels of the image divided by ES. Unchanged 
on exit. 


WE. SIZE 
On entry WE_SIZE contains the width of the image divided by ES. Unchanged on exit. 


IMAGE.SUM 
On exit IMAGE_SUM is the sum of IMAGE and SCAL modulo 256. 


Errors 


No explicit error checking is done. 


Auxiliary Routines 


None. 


Accuracy 
Not applicable. 
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9 Further Comments 


None. 


10 Keywords 


Image Processing 
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4.3 SUB_8 
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Purpose 


This routine subtracts two 8-bit two’s complement images. The result is returned modulo 
206. 


Specification 
SUB. 8 (IMAGE1, IMAGE2,NS_SIZE,WE_SIZE,IMAGE- DIFF) 


INTEGER*1 IMAGE1 (,, NS_SIZE,WE-SIZE), 
& IMAGE2 (, ,NS_SIZE,WE-_SIZE), 
& IMAGE. DIFF (, , NS_SIZE,WE- SIZE) 


INTEGER*4 NS_SIZE,WE_SIZE 


Description 
IMAGE2 is subtracted modulo 256 from IMAGE1 pixel-by-pixel to form IMAGE. DIFF. 


References 
None. 


Arguments 
IMAGEI1 


On entry is one of the 8-bit images to be subtracted. Each dimension of the image 
size must be >HSand a multiple of ES. IMAGE1 may be overwritten by the result if 
desired. 


IMAGE2 


On entry IMAGE2 is the other image to be subtracted. The dimensions of IMAGE2 
must be identical to those of IMAGE1. IMAGE2 may be overwritten by the result if 
desired. 


NS.SIZE 


On entry NS_SIZE contains the height in pixels of the images divided by ES. Unchanged 
on exit. 


WE. SIZE 


On entry WE-.SIZE contains the width of the images divided by #.S. Unchanged on 
exit. 


IMAGE. DIFF 
On exit IMAGE_DIFF is IMAGE1 minus IMAGE2 modulo 256. 


Errors 


No explicit error checking is done. 


Auxiliary Routines 
None. 
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8 Accuracy 
Not applicable. 


9 Further Comments 


None. 


10 Keywords 


Image Processing 


Image Processing Library 
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MULT. 8. TO_ 16 Image Processing Primitives 





4.4 MULT_8_TO_16 
1 Purpose 
This routine multiplies two 8-bit two’s complement images to form a 16-bit image. 
2 Specification 
MULT_8.TO-_16 (IMAGE1,IMAGE2,NS_SIZE,WE-.SIZE,IMAGE. PROD) 


INTEGER*1 IMAGE1 (, ,NS_SIZE,WE- SIZE), 
& IMAGE2 (, ,NS_SIZE,WE- SIZE) 


INTEGER*2 IMAGE. PROD (, , NS_SIZE,WE- SIZE) 
INTEGER*4 NS_SIZE,WE-_SIZE 
3 Description 
IMAGE1 is multiplied by IMAGE2 pixel-by-pixel to form IMAGE_ PROD. 


4 References 


None. 





5 Arguments 
IMAGE1 


On entry is one of the 8-bit images to be multiplied together. Each dimension of the 
image size must be > Sand a multiple of ES. Unchanged on exit. 


IMAGE2 


On entry IMAGE? is the other image to be multiplied. The dimensions of IMAGE2 must 
be identical to those of IMAGE1. Unchanged on exit. 


NS_SIZE 


On entry NS_SIZE contains the height in pixels of the images divided by ES . Unchanged 
on exit. 


WE-.SIZE 


On entry WE_SIZE contains the width of the images divided by ES. Unchanged on 
exit. 


IMAGE. PROD 
On exit IMAGE. PROD is the product of IMAGE1 and IMAGE2. 


6 Errors 


No explicit error checking is done. 


7 Auxiliary Routines 


None. 
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8 Accuracy 
Not applicable. 


9 Further Comments 


None. 


10 Keywords 


Image Processing 
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4.5 SCAL.MULT_8_TO_16 


42 


Purpose 


This routine multiplies an 8-bit two’s complement image by an 8-bit two’s complement scalar 
to form a 16-bit image. 


Specification 
SCAL. MULT_8_TO_16 (IMAGE,SCAL,NS_SIZE,WE-_SIZE,IMAGE. PROD) 


INTEGER*1 IMAGE (,, NS_SIZE,WE_SIZE), 
& SCAL 


INTEGER*2 IMAGE. PROD (,, NS_SIZE,WE. SIZE) 


INTEGER*4 NS_SIZE,WE_SIZE 


Description 
IMAGE is multiplied by SCAL pixel-by-pixel to form IMAGE. PROD. 


References 


None. 


Arguments 


IMAGE 


On entry IMAGE is the 8-bit image to be multiplied by SCAL. Each dimension of the 
image size must be >HSand a multiple of ES. Unchanged on exit. 


SCAL 
On entry SCAL is the scalar that multiplies IMAGE. Unchanged on exit. 


NS.SIZE 


On entry NS_SIZE contains the height in pixels of the image divided by ES. Unchanged 
on exit. 


WE. SIZE 


On entry WE.SIZE contains the width of the image in pixels divided by ES. Unchanged 
on exit. 


IMAGE. PROD 
On exit IMAGE. PROD is the product of IMAGE and SCAL. 


Errors 


No explicit error checking is done. 


Auxiliary Routines 


None. 
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8 Accuracy | 
Not applicable. 


9 Further Comments 


None. 


10 Keywords 


Image Processing 
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4.6 SCAL_DIV_8 
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it 


Purpose 


This routine divides an 8-bit two’s complement image by an 8-bit two’s complement scalar 
to form an 8-bit image. 


Specification 
SCAL_DIV_8 (IMAGE,SCAL,NS_SIZE,WE.SIZE, IMAGE. QUOT) 
INTEGER*1 IMAGE (,, NS_SIZE,WE_SIZE), 
& SCAL, | 
& IMAGE. QUOT (, , NS_SIZE,WE-_SIZE) 
INTEGER*4 NS_SIZE,WE-SIZE 
Description 


IMAGE is divided by SCAL pixel-by-pixel to form IMAGE. QUOT. 


References 


None. 


Arguments 


IMAGE 


On entry IMAGE is the 8-bit image to be divided by SCAL. Each dimension of the image 
size must be >ESand a multiple of ES. Unchanged on exit. 


SCAL 
On entry SCAL is the scalar that divides IMAGE. Unchanged on exit. 


NS.SIZE 


On entry NS_SIZE contains the height in pixels of the image divided by ES. Unchanged 
on exit. 


WE. SIZE 
On entry WE-SIZE contains the width of the image divided by ES. Unchanged on exit. 


IMAGE. QUOT 
On exit IMAGE_ QUOT is the quotient of IMAGE and SCAL. 


Errors 


No explicit error checking is done. 


Auxiliary Routines 


None. 
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8 Accuracy 
Not applicable. 


9 Further Comments 
None. 


10 Keywords 


Image Processing 
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Shifts 


4.7 SHIFT_IMAGE_NORTH_P 
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Purpose 


SHIFT_IMAGE.NORTH-P shifts an entire image north with planar boundary conditions. 
Zeroes are introduced at the boundaries. 


Specification 


SUBROUTINE SHIFT_IMAGE_NORTH-P (IMAGE, BITS, NS_SIZE, WE-SIZE, 
& SHIFT. IMAGE, DIST) 


LOGICAL IMAGE (,, BITS,NS.SIZE,WE-_SIZE), 
& SHIFT.IMAGE (,, BITS,NS.SIZE,WE_SIZE) 7 


INTEGER*4 NS_SIZE, WE.SIZE, DIST, BITS 


Description 


All of IMAGE is shifted north by DIST pixels with planar boundary conditions. The precision 
of IMAGE is specified in BITS. 


References 
None. 


Arguments 
IMAGE 


On entry IMAGE is the image to be shifted. Each dimension of the image size must be 
in the range ESto 32ESand a multiple of ES. The image is assumed to be stored in 
IMAGE using a crinkled mapping where the third and fourth co-ordinates of IMAGE 
are vertical and horizontal co-ordinates of each crinkled tile. Unchanged on exit, unless 
IMAGE and SHIFT.IMAGE occupy the same area of store. 


BITS 
On entry BITS is the precision in bits of the image. Unchanged on exit. 


NS-SIZE 


On entry NS.SIZE contains the height of IMAGE in pixels divided by ES. NS_SIZE © 
must be in the range 1 to 32. Unchanged on exit. 


WE. SIZE 


On entry WE-SIZE contains the width of IMAGE in pixels divided by ES. WE_SIZE 
must be in the range 1 to 32. Unchanged on exit. 
SHIFT. IMAGE 


On exit SHIFT. IMAGE contains the shifted version of IMAGE. SHIFT.IMAGE may 
occupy the same area of store as IMAGE. (They must occupy identical or disjoint but 
not overlapping areas of store.) 
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DIST 


On entry DIST is the distance IMAGE is to be moved. DIST must be non-negative. 
Unchanged on exit. 


6 Errors 


No explicit error checking is done. 


7 Auxiliary Routines 


None. 


8 Accuracy 
Not applicable. 


9 Further Comments 


260 stack planes are allocated. BITS x DISTmoduloNS_SIZE planes will be needed. 
For east or west shifts BITS x DISTmoduloW E_SIZE planes will be needed. For large 
shifts to be performed you will need to increase the program’s stack allocation (see stack size 
in DAP Series: Program Development under UNIX or DAP Series: Program Development 
under VAX/VMS). . 


10 Keywords 


Shifting, Image Processing 
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4.8 SHIFT_IMAGE_NORTH_C 
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Purpose 


SHIFT.IMAGE_NORTH_C shifts an entire image north with cyclic (wrap-around) bound- 
ary conditions. 


Specification 


SUBROUTINE SHIFT_IMAGE_NORTH_C (IMAGE, BITS, NS_SIZE, WE-SIZE, 
& SHIFT_IMAGE, DIST) 


LOGICAE IMAGE (,, BITS,NS_SIZE,WE-_SIZE), 
& SHIFT_IMAGE (,, BITS,NS_SIZE,WE-SIZE) 


INTEGER*4 NS_SIZE, WE_SIZE, DIST, BITS 


Description 


All of IMAGE is shifted north by DIST pixels with cyclic boundary conditions. The precision 
of IMAGE is specified in BITS. 


References 


None. 


Arguments 


IMAGE 


On entry IMAGE is the image to be shifted. Each dimension of the image size must be 
in the range ESto 32ES and a multiple of ES. The image is assumed to be stored in 
IMAGE using a crinkled mapping where the third and fourth co-ordinates of IMAGE 
are vertical and horizontal co-ordinates of each crinkled tile. Unchanged on exit, unless 
IMAGE and SHIFT_IMAGE occupy the same area of store. 


BITS 
On entry BITS is the precision in bits of the image. Unchanged on exit. 


NS_SIZE 


On entry NS_SIZE contains the height in pixels of IMAGE divided by ES. NS_SIZE 
must be in the range 1 to 32. Unchanged on exit. 


WE. SIZE 


On entry WE-_SIZE contains the width in pixels of IMAGE divided by ES. WE_SIZE 
must be in the range 1 to 32. Unchanged on exit. 


SHIFT.IMAGE 


On exit SHIFT_IMAGE contains the shifted version of IMAGE. SHIFT.IMAGE may 
occupy the same area of store as IMAGE. (They must occupy identical or disjoint but 
not overlapping areas of store.) 


DIST 


On entry DIST is the distance IMAGE is to be moved. DIST must be non-negative. 
Unchanged on exit. 
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6 Errors 


No explicit error checking is done. 


7 Auxiliary Routines 


None. 


8 Accuracy 
Not applicable. 


9 Further Comments 


260 stack planes are allocated. BITS x DISTmoduloNS_SIZE planes will be needed. 
For east or west shifts BITS x DISTmoduloW E_SIZE planes will be needed. For large 
shifts to be performed you will need to increase the program’s stack allocation (see stack size 
in DAP Series: Program Development under UNIX or DAP Series: Program Development 
under VAX/VMS). 


10 Keywords 


Shifting, Image Processing 
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4.9 SHIFT.ROW_NORTH_P 
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Purpose 


SHIFT: ROW. NORTH_P shifts a row of matrices of an image north with planar boundary 
conditions. Zeroes are introduced at the boundaries. 


Specification 


SUBROUTINE SHIFT. ROW.NORTH_P (IMAGE, BITS, NS_SIZE, WE-SIZE, 
& SHIFT. ROW, ROW, DIST) 


LOGICAL IMAGE (,, BITS,NS._SIZE,WE-_SIZE), 
& SHIFT. ROW (,, BITS,WE_SIZE) 


INTEGER*4 NS_SIZE, WE_SIZE, ROW, DIST, BITS 


Description 


SHIFT. ROW is the ROW* row of DAP matrices from the version of IMAGE shifted north 
by DIST pixels with planar boundary conditions. Thus if IMAGE1 were the (imaginary) 
shifted version of IMAGE, then SHIFT. ROW (, ,J)=IMAGE1(,, ROW, J) for J=1,2,...WE.SIZE. 


References 


None. 


Arguments 


IMAGE 


On entry IMAGE is the image to be shifted. Each dimension of the image size must be 
in the range ESto 32ES and a multiple of ES. The image is assumed to be stored in 
IMAGE using a crinkled mapping where the third and fourth co-ordinates of IMAGE 
are vertical and horizon- tal co-ordinates of each crinkled tile. Unchanged on exit. 


BITS 


On entry BITS is the precision in bits of the image. Unchanged on exit. 


NS.SIZE 


On entry NS.SIZE contains the height in pixels of IMAGE divided by ES. NS_SIZE 
must be in the range 1 to 32. Unchanged on exit. 


WE-.SIZE 


On entry WE_SIZE contains the width in pixels of IMAGE divided by ES. WE_SIZE 
must be in the range 1 to 32. Unchanged on exit. 
SHIFT. ROW 


On exit SHIFT. ROW contains the ROW“ row of DAP matrices of the shifted version 
of IMAGE. 


ROW 
ROW is the row of the shifted image to be calculated. 
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DIST 


On entry DIST is the distance IMAGE is to be moved. DIST must be non-negative. 
Unchanged on exit. 


Errors 


No explicit error checking is done. 


Auxiliary Routines 


None. 


Accuracy 
Not applicable. 


Further Comments 


260 stack planes are allocated. BITS x DISTmoduloNS_SIZE planes will be needed. 
For east or west shifts BITS x DISTmoduloW E_SIZE planes will be needed. For large 
shifts to be performed you will need to increase the program’s stack allocation (see stack size 
in DAP Series: Program Development under UNIX or DAP Series: Program Development 
under VAX/VMS). 


Keywords 


Shifting, Image Processing 


4.10 SHIFT.ROW_NORTH_C 
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Purpose 


SHIFT.ROW_NORTH.C shifts a row of matrices of an image north with cyclic (wrap- 
around) boundary conditions. 


Specification 


SUBROUTINE SHIFT. ROW_NORTH_C (IMAGE, BITS, NS_SIZE, WE-SIZE, 
& SHIFT. ROW, ROW, DIST) 


LOGICAL IMAGE(, , BITS,NS_SIZE,WE.SIZE), SHIFT- ROW (,, BITS, WE_SIZE) 
INTEGER*4 NS_SIZE, WE.SIZE, ROW, DIST, BITS 


Description 


SHIFT. ROW is the ROW** row of DAP matrices from the version of IMAGE shifted north 
by DIST pixels with cyclic boundary conditions. Thus if IMAGE1 were the (imaginary ) 
shifted version of IMAGE, then SHIFT. ROW (, ,J)=IMAGE1 (,, ROW, J) for J=1,2,...WE_SIZE. 


References 
None. 
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5 Arguments 
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IMAGE 


On entry IMAGE is the image to be shifted. Each dimension of the image size must be 
in the range HS to 32HS and a multiple of ES. The image is assumed to be stored in 
IMAGE using a crinkled mapping where the third and fourth co-ordinates of IMAGE 
are vertical and horizontal co-ordinates of each crinkled tile. Unchanged on exit. 


BITS 


On entry BITS is the precision in bits of the image. Unchanged on exit. . 


NS_SIZE 


On entry NS_SIZE contains the height in pixels of IMAGE divided by ES. NS_SIZE 
must be in the range 1 to 32. Unchanged on exit. 


WE. SIZE 


On entry WE_SIZE contains the width in pixels of IMAGE divided by ES. WE_SIZE 
must be in the range 1 to 32. Unchanged on exit. 


SHIFT.ROW 


On exit SHIFT. ROW contains the ROW“ row of DAP matrices of the shifted version 
of IMAGE. 
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ROW 
ROW is the row of the shifted image to be calculated. 
DIST 


On entry DIST is the distance IMAGE is to be moved. DIST must be non-negative. 
Unchanged on exit. 


6 Errors 


No explicit error checking is done. 


7 Auxiliary Routines 


None. 


8 Accuracy 
Not applicable. 


9 Further Comments 


260 stack planes are allocated. BITS x DISTmoduloNS_SIZE planes will be needed. 
For east or west shifts BITS x DISTmoduloW E_SIZE planes will be needed. For large 
shifts to be performed you will need to increase the program’s stack allocation (see stack size 
in DAP Series: Program Development under UNIX or DAP Series: Program Development 
under VAX/VMS). 





10 Keywords | 
Shifting, Image Processing 
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4.11 SHIFT_.COL_.NORTH_P 


1 Purpose 


SHIFT. COL. NORTH-P shifts a column of matrices of an image north with planar boundary 
conditions. Zeroes are introduced at the boundaries. 


2 Specification 


SUBROUTINE SHIFT_COL.NORTH-_P (IMAGE, BITS, NS_SIZE, WE-SIZE, 
& SHIFT_COL, COL, DIST) 


LOGICAL IMAGE(, , BITS,NS.SIZE,WE_SIZE), SHIFT_COL (, , BITS,NS_SIZE) 
INTEGER*4 NS_SIZE, WE_SIZE, COL, DIST, BITS 


3 Description 


SHIFT_COL is the COL column of DAP matrices from IMAGE shifted north by DIST 
pixels with planar boundary conditions. Thus if IMAGE1 were the (imaginary) shifted 
version of IMAGE, then SHIFT.COL(,, I=IMAGE1(, ,I,COL) for I=1, 2,...NS_SIZE. 


4 References 


None. 





5 Arguments 
IMAGE 


On entry IMAGE is the image to be shifted. Each dimension of the image size must be 
in the range ESto 32S and a multiple of ES. The image is assumed to be stored in 
IMAGE using a crinkled mapping where the third and fourth co-ordinates of IMAGE 
are vertical and horizontal co-ordinates of each crinkled tile. Unchanged on exit. 


BITS 


On entry BITS is the precision in bits of the image. Unchanged on exit. | 


NS.SIZE 


On entry NS_SIZE contains the height in pixels of IMAGE divided by ES. NS_SIZE 
must be in the range 1 to 32. Unchanged on exit. 


WE. SIZE 


On entry WE-SIZE contains the width in pixels of IMAGE divided by ES. WE_SIZE 
must be in the range 1 to 32. Unchanged on exit. 


SHIFT.COL 


On exit SHIFT. COL contains the COL” column of DAP matrices of the shifted version 
of IMAGE. 


COL 
COL is the column of the shifted image to be calculated. 
DIST 


On entry DIST is the distance IMAGE is to be moved. DIST must be non-negative. 
Unchanged on exit. 
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Errors 


No explicit error checking is done. 


Auxiliary Routines 
None. 


Not available. 


Accuracy 
Not applicable. 


Further Comments 


260 stack planes are allocated. BITS x DISTmoduloNS.SIZE planes will be needed. 
For east or west shifts BITS x DISTmoduloW E_SIZE planes will be needed. For large 
shifts to be performed you will need to increase the program’s stack allocation (see stack size 
in DAP Series: Program Development under UNIX or DAP Series: Program Development 
under VAX/VMS). 


Keywords 
Shifting, Image Processing 
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4.12 SHIFT_COL.NORTH_C 


1 Purpose 


SHIFT.COL.NORTH.C shifts a column of matrices of an image north with cyclic (wrap- 
around) boundary conditions. 


2 Specification 


SUBROUTINE SHIFT_COL.NORTH_C (IMAGE, BITS, NS_SIZE, WE_SIZE, 
& SHIFT. COL, COL, DIST) 


LOGICAL IMAGE(, , BITS,NS_SIZE,WE_SIZE), SHIFT. COL(, , BITS,NS_SIZE) 
INTEGER*4 NS_SIZE, WE.SIZE, COL, DIST, BITS 


3 Description 


SHIFT.COL is the COL** column of DAP matrices from IMAGE shifted north by DIST 
pixels with cyclic boundary conditions. Thus if IMAGE1 were the (imaginary) shifted version 
of IMAGE, then SHIFT.COL(,,ID=IMAGE1(,,1I,COL) for I=1,2,...NS_SIZE. 


4 References 
None. 





5 Arguments 
IMAGE 


On entry IMAGE is the image to be shifted. Each dimension of the image size must 
be > ESand a multiple of ES. The image is assumed to be stored in IMAGE using 
a crinkled mapping where the third and fourth co-ordinates of IMAGE are vertical and 
horizontal co-ordinates of each crinkled tile. Unchanged on exit. 


BITS 
On entry BITS is the precision in bits of the image. Unchanged on exit. 


NS.SIZE 


On entry NS_SIZE contains the height in pixels of IMAGE divided by ES. NS_SIZE 
must be in the range 1 to 32. Unchanged on exit. 


WE-SIZE 


On entry WE_SIZE contains the width in pixels of IMAGE divided by ES. WE_SIZE 
must be in the range 1 to 32. Unchanged on exit. 


SHIFT.COL 


On exit SHIFT. COL contains the COL** column of DAP matrices of the shifted version 
of IMAGE. 


COL 
COL is the column of the shifted image to be calculated. 
DIST 


On entry DIST is the distance IMAGE is to be moved. DIST must be non-negative. 
Unchanged on exit. 
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Errors 


No explicit error checking is done. 


Auxiliary Routines 
None. 


Accuracy 
Not applicable. 


Further Comments 


260 stack planes are allocated. BITS x DISTmoduloNS_SIZE planes will be needed. 
For east or west shifts BITS x DISTmoduloW E_.SIZE planes will be needed. For large 
shifts to be performed you will need to increase the program’s stack allocation (see stack size 


in DAP Series: Program Development under UNIX or DAP Series: Program Development 
under VAX/VMS). 


Keywords 


Shifting, Image Processing 


Image Processing Library man014.03 57 


SHIFT.SHEET. NORTH. P Image Processing Primitives 





4.13 SHIFT_SHEET_NORTH_P 


1 Purpose 


SHIFT.SHEET_.NORTH_P shifts a matrix of an image north with planar boundary condi- 
tions. Zeroes are introduced at the boundaries. 


2 Specification 


SUBROUTINE SHIFT.SHEET. NORTH. P (IMAGE, BITS, NS_SIZE, WE-_SIZE, 
& SHIFT_SHEET, ROW, COL, DIST) 


LOGICAL IMAGE (,, BITS,NS_SIZE,WE_SIZE), SHIFT. SHEET (, , BITS) 


INTEGER*4 NS_SIZE, WE_SIZE, ROW, COL, DIST, BITS 


3 Description 


SHIFT.SHEET is the (ROW ,COL**) DAP matrix from IMAGE if it were shifted north 
by DIST pixels with planar boundary conditions. Thus if IMAGE1 were (the imaginary) 
shifted version of IMAGE, then SHIFT.SHEET (, )=IMAGE1(,, ROW, COL). 


4 References 


None. 





5 Arguments 
IMAGE 


On entry IMAGE is the image to be shifted. Each dimension of the image size must be 
in the range ESto 32ES and a multiple of ES. The image is assumed to be stored in 
IMAGE using a crinkled mapping where the third and fourth co-ordinates of IMAGE 
are vertical and horizontal co-ordinates of each crinkled tile. Unchanged on exit. 


BITS 


On entry BITS is the precision in bits of the image. Unchanged on exit. 


NS.SIZE 


On entry NS_SIZE contains the height in pixels of IMAGE divided by ES. NS_SIZE 
must be in the range 1 to 32. Unchanged on exit. 


WE. SIZE 


On entry WE_SIZE contains the width in pixels of IMAGE divided by ES. WE_SIZE 
must be in the range 1 to 32. Unchanged on exit. 


SHIFT.SHEET 


On exit SHIFT.SHEET contains the (ROW** ,COL) DAP matrix of the shifted version 
of IMAGE. 
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ROW 
ROW is the row of the shifted image to be calculated. 


DIST 


On entry DIST is the distance IMAGE is to be moved. DIST must be non-negative. 
Unchanged on exit. 


6 Errors 


No explicit error checking is done. 


7 Auxiliary Routines 


None. 


8 Accuracy 
Not applicable. 


9 Further Comments 


None. 


10 Keywords 


Shifting, Image Processing 
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4.14 SHIFT_SHEET_NORTH_C 


1 Purpose 


SHIFT. SHEET. NORTH-C shifts a matrix of an image north with cyclic (wrap-around) 
boundary conditions. 


2 Specification 


SUBROUTINE SHIFT_SHEET. NORTH_C (IMAGE, BITS, NS_SIZE, WE-_SIZE, 
& SHIFT.SHEET, ROW, COL, DIST) 


LOGICAL IMAGE(, , BITS,NS_SIZE,WE_SIZE), SHIFT. SHEET (,, BITS) 
INTEGER*4 NS_SIZE, WE_SIZE, ROW, COL, DIST, BITS 


3 Description | 


SHIFT.SHEET is the (ROW** ,COL**) DAP matrix from IMAGE if it were shifted north by 
DIST pixels with cyclic boundary conditions. Thus if IMAGE1 were the (imaginary) shifted 
version of IMAGE, then SHIFT.SHEET ( ,)=IMAGE1(,, ROW, COL). 


4 References 


None. 





5 Arguments 
IMAGE 


On entry IMAGE is the image to be shifted. Each dimension of the image size must be 
in the range ESto 32ESand a multiple of ES. The image is assumed to be stored in 
IMAGE using a crinkled mapping where the third and fourth co-ordinates of IMAGE 
are vertical and horizontal co-ordinates of each crinkled tile. Unchanged on exit. 


BITS 
On entry BITS is the precision in bits of the image. Unchanged on exit. 


NS.SIZE 


On entry NS_SIZE contains the height in pixels of IMAGE divided by ES. NS_SIZE 
must be in the range 1 to 32. Unchanged on exit. 


WE-_SIZE 


On entry WE_SIZE contains the width in pixels of IMAGE divided by ES. WE_SIZE 
must be in the range 1 to 32. Unchanged on exit. 


SHIFT.SHEET 


On exit SHIFT_SHEET contains the (ROW ,COL'*) DAP matrix of the shifted version 
of IMAGE. 


ROW 
ROW is the row of the shifted image to be calculated. 
DIST 


On entry DIST is the distance IMAGE is to be moved. DIST must be non-negative. 
Unchanged on exit. 
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6 Errors 


No explicit error checking is done. 


7 Auxiliary Routines 


None. 


8 Accuracy 
Not applicable. 


9 Further Comments 


None. 


10 Keywords 


Shifting, Image Processing 
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Chapter 5 


Low Level Image Processing 
Routines 


These routines perform low level image processing tasks such as convolutions and FFT’s. Routines 
specified as 8-bits per pixel may be used for 16-bits per pixel images by replacing the 8 in the name 
by 16. For example: 


AVERAGE_8 — convolves a square averaging mask with an 8-bit two’s complement 
image. 

AVERAGE. 16 —- convolves a square averaging mask with a 16-bit two’s complement 
image. 


The COMP._FFT_2D routines which perform FFTs on images with 8, 16 and 24 bits per pixel, 
are listed in full. 
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5.1 ABS. THRESH_8 


1 Purpose 
ABS. THRESH. 8 thresholds an 8-bit two’s complement image. 


2 Specification 


SUBROUTINE ABS. THRESH_8 (IMAGE , NS_SIZE, WE_SIZE, THRESH, 
& HIGH.TIDE, FLIP, IFAIL) 


INTEGER*1 IMAGE(, , NS_SIZE,WE-SIZE) 
INTEGER*4 NS_SIZE, WE_SIZE, HIGH_TIDE, IFAIL 


LOGICAL THRESH (, , NS.SIZE, WE_SIZE) , FLIP 


3 Description 


THRESH is a logical bit-map of the image which has .TRUE. values where the IMAGE 
grey-scale values exceed HIGH. TIDE and .FALSE. otherwise. (This is reversed if FLIP is 
-FALSE. instead of .TRUE.) 


4 References 


None. 


5 Arguments 
IMAGE 


On entry IMAGE contains the image to be thresholded. Each dimension of the image 
size must be in the range Sto 32HSin multiples of ES. The image is assumed to be 
stored in IMAGE using a crinkled mapping where the third and fourth co-ordinates of 
IMAGE are the vertical and horizontal co-ordinates of each crinkled tile. Unchanged on 
exit. 

NS_SIZE 


On entry NS_SIZE contains the height in pixels of the image divided by FS. NS_SIZE 
must be in the range 1-32. Unchanged on exit. 


WE. SIZE 


On entry WE_SIZE contains the width of the image in pixels divided by ES. WE_SIZE 
must be in the range 1-32. Unchanged on exit. 


THRESH 


On exit THRESH contains a logical bit-map of from the thresholding of IMAGE. 
THRESH stores the bit-map using a crinkled mapping (see IMAGE above). 


HIGH. TIDE 
HIGH-TIDE is the threshold with which IMAGE is compared. 
FLIP 


If FLIP is .TRUE., THRESH is .TRUE. when IMAGE is greater than HIGH. TIDE. If 
FLIP is .FALSE., THRESH is .FALSE. when IMAGE is greater than HIGH_TIDE. 
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IFAIL 


IF AIL has the value 0 on exit unless an error occurs. (See 6 below.) 


6 Errors 
IFAIL=0 Successful exit 
TFAIL=1 NS_SIZE out of range 
IFAIL=2 WE_SIZE out of range 
IFAIL=13 HIGH. TIDE out of range —128 to 127 


7 Auxiliary Routines 


None. 


8 Accuracy 
Not applicable. 


9 Further Comments 


None. 


10 Keywords 


Thresholding, Image Processing 
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AVERAGE_8 Low Level Image Processing Routines 


5.2 AVERAGE_8 


1 Purpose 


AVERAGE.8 convolves a square averaging mask with an 8-bit two’s complement image. 


2 Specification 


SUBROUTINE AVERAGE. 8 (IMAGE,NS_SIZE, WE_SIZE, AVG, 
& MASK _SIZE, IFAIL) 


INTEGER*1 IMAGE(,,NS_SIZE,WE_SIZE), AVG(, , NS_SIZE,WE-SIZE) 


INTEGER*4 NS_SIZE, WE.SIZE, MASK_SIZE, IFAIL 


3 Description 


AVG is the convolution of IMAGE and a centred square averaging mask. IMAGE is assumed 
to be extended by a field of zeros for the convolution. 


4 References. 


None. 


5 Arguments 
IMAGE 


On entry IMAGE contains the image to be convolved with the averaging mask. Each 
dimension of the image size must be in the range ESto 32ESin multiples of HS. The 
image is assumed to be stored in IMAGE using a crinkled mapping where the third 
and fourth co-ordinates of IMAGE are the vertical and horizontal co-ordinates of each 
crinkled tile. 


IMAGE is overwritten. 
NS.SIZE 


On entry NS_SIZE contains the height in pixels of the image divided by ES. NS_SIZE 
must be in the range 1-32. Unchanged on exit. 


WE._SIZE 


On entry WE_SIZE contains the width of the image divided by HS. WE-SIZE must be 
in the range 1-32. Unchanged on exit. 


AVG 


On exit AVG contains the convolution of IMAGE with an averaging mask. AVG stores 
the convolution using a crinkled mapping (see IMAGE above). 


MASK_SIZE 


On entry MASK.SIZE must be one of 3, 5, 7,...,31. This is length of the side of the 
averaging mask. Unchanged on exit. 


IFAIL 


IFAIL has the value 0 on exit unless an error occurs. (See 6 below.) 
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6 Errors 
IFAIL=0 Successful exit 
IFAIL=1 NS.SIZE out of range 
IFAIL=2 WE.SIZE out of range 
IFAIL=3 . MASK._SIZE invalid 


7 Auxiliary Routines 


Image Processing Library routines SHIFT.COL. WEST-P, SHIFT.COL_EAST-P, 
SHIFT. ROW. WEST. P and SHIFT. ROW_EAST-_P are called. 


8 Accuracy 
Not available. 


9 Further Comments 


None. : 


10 Keywords 


Convolution, Local Averaging, Image Processing 








Image Processing Library man014.03 67 





BOX_IN. BOX. 8 Low Level Image Processing Routines 





5.3 BOX_IN_.BOX_8 


1 Purpose 


BOX_IN. BOX.8 applies box-in-a-box prescreening to an 8-bit per pixel two’s complement 
image. The result is compared with a user supplied threshold and a binary result returned. 


2 Specification 


SUBROUTINE BOX.IN. BOX.8 (IMAGE, NS_SIZE, WE.SIZE, BOX, 


& BOX_SIZE, THRESH, W1, W2, IFAIL) 
INTEGER*1 IMAGE(,,NS.SIZE,WE_SIZE), 

& W1(,,NS_SIZE,WE-SIZE), 

& W2(,, NS_SIZE,WE_SIZE) 


LOGICAL BOX(,,NS.SIZE, WE.SIZE) 


INTEGER*4 NS_SIZE, WE.SIZE, BOX.SIZE, THRESH, IFAIL 


3 Description 


For each pixel in IMAGE the grey-scale values in a square box centred on the pixel are 
averaged. An average of the grey-scale values is also calculated for a box of approximately 
one-half the area of the first box. The average over the large box is subtracted from the 
average over the small box and the result is compared with THRESH. If the result is greater 
than THRESH, BOX will be .TRUE otherwise BOX will be .FALSE.. 





4 References 


None. 


5 Arguments 


IMAGE 


On entry IMAGE contains the image to be pre-screened. Each dimension of IMAGE 
must be in the range ES to 32ES in multiples of HS . The image is assumed to be stored 
using a crinkled mapping where the third and fourth co-ordinates of IMAGE are the 
vertical and horizontal co-ordinates of each crinkled tile. 


IMAGE is overwritten. 
NS.SIZE 


On entry NS_SIZE contains the height in pixels of the image divided by ES. NS_SIZE 
must be in the range 1 to 32. Unchanged on exit. 


WE. SIZE 


On entry WE_SIZE contains the width of the image in pixels divided by ES. WE_SIZE 
must be in the range 1 to 32. Unchanged on exit. 





68 man014.03 AMT 











Low Level Image Processing Routines BOX. IN_ BOX. 8 


10 


BOX 


On exit BOX contains the result the thresholding operation stored using a crinkled 
mapping. 


BOX_SIZE 


On entry BOX_SIZE is the size of the square outer box in pixels. It must be one of 5, 7, 
9,...31. Unchanged on exit. 


THRESH 


On entry THRESH is the threshold with which the difference of the two averages is 
compared. Unchanged on exit. 


Wi 
Workspace. 
W2 
Workspace. 
IFAIL 


IFAIL equals 0 on exit unless an error occurs. (See 6 below) 


Errors 
IFAIL=0 Successful exit 
IFAIL=1 NS-.SIZE out of range 
IFAIL=2 WE.SIZE out of range 
IFAIL=3 © BOX.SIZE out of range 
Auxiliary Routines 


The Image Processing Library routines AVERAGE_8, SHIFT.COL.WEST-_P, 
SHIFT. COL_EAST-_P, SHIFT. ROW. WEST_P and SHIFT. ROW_EAST.P are called. 


Accuracy 


Not available. 


Further Comments 


None. 


Keywords 


Prescreening, Image Processing 
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5.4 C0O6_FFT_ESS 


1 Purpose 


C06_FFT_ ESS calculates the two dimensional discrete Fourier transform of HS x ES complex 
points. 


2 Specification 
SUBROUTINE C06_FFT_ESS( X , Y , INVERS , FIRST ) 
REAL X(,) , Y(,) 


LOGICAL INVERS , FIRST 


3 Description 


The two dimensional transform is calculated by performing independent sets of row and 
column ES-point transforms. 


The data is then in bit reversed order independently in rows and columns, and a final shuffle 
is performed to reorder the data. 


For a description of the general theory of FFTs see [1]. 





4 References 


(1] BRIGHAM E.O. 
The Fast Fourier Transform. 
Prentice-Hall, 1974 


5 Arguments 
X —~ REAL MATRIX 


On entry X contains the real part of the data to be transformed. On exit X contains the 
real part of the transformed data. 


Y - REAL MATRIX 


On entry Y contains the imaginary part of the data to be transformed. On exit Y 
contains the imaginary part of the transformed data. 


INVERS — LOGICAL 
If INVERS is set to .FALSE. the transform: 


Xje + tYjR - > > (Amn + {Bran eaponi— im — 7) <a (& = 1)(n = 1) 


ES 
ms n 
is calculated, where j = 1,2,... , BS;k =1,2,... , ES and the summations are also 
overm=1,2,..., HSandn=1,2,... , ES. 





70 man014.03 AMT 














Low Level Image Processing Routines C06_FFT_ESS 


If INVERS is set to .TRUE. the transform: 


: (m—l1)(j7-1 —1)(k ~1 
Annie. = (Xj + Yjx)eap-2ri 2 — WG —Y) + ene a) 
jk 


is calculated, where m = 1,2,... , HS;n=1,2,... , ES and the summations are also 
over j= 1,2,..., HSand k= 1,2,... , ES 


FIRST — LOGICAL 


If FIRST is set to .TRUE. the exponential coefficients for the transform are calculated. 
Consequently FIRST must be set to .TRUE. the first time this routine is called within 
a program, but may be set to .FALSE. for all subsequent calls. 


6 Errors 


None. 


7 Auxiliary Routines 


This routine calls the DAP library routines Z.C06_F2DCOEFF, Z_CO6_ROWFFT, 
Z-C06_COLFFT and Z.C06_F2DBREV. 


8 Accuracy 


Accuracy will be data dependent. Some indication of the accuracy may be obtained by 
performing a subsequent inverse transform and comparing the results with the original data. 


9 Further Comments 


This routine uses a common block with the name CCO6FFTESSQ. Consequently the user 
program must not use a common block with this name. 


10 Keywords 


Fast Fourier Transform 


11 Example 


Not available. 
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5.56 CO6_FFT_LV 


1 Purpose 
CO6_FFT_LV performs a one dimensional finite Fourier transform of HS? complex points. 


2 Specification 


SUBROUTINE CO6_FFT_LV( X,Y , INVERS , FIRST ) 
REAL X(,) , Y() 
LOGICAL INVERS , FIRST 


3 Description 


The data is considered as HS? complex points in long vector order, and the transform is 
calculated by performing linked row and column transforms. The first step is to calculate 
ES-point transforms along each row of complex data. The results of the row transforms 
are multiplied by a second set of exponential factors and then HS-point transforms are 
calculated along each column in a similar way to the row transforms but using different 
exponential factors. The exponential factors are set up in such a way as to ensure that 
the row and column transforms are linked correctly to give the required one dimensional 
‘transform. The final step reorders the data which is in bit reversed order. 


For a description of the general theory of FFTs see [1]. 





4 References 


) [1] BRIGHAM E.O. 
| The Fast Fourier Transform. 
| Prentice-Hall, 1974 


5 Arguments 
X — REAL MATRIX 


On entry X contains the real part of the data to be transformed. On exit X contains the 
transformed real part of the data. 


Y —- REAL MATRIX 


On entry Y contains the imaginary part of the data to be transformed. On exit Y 
contains the transformed imaginary part of the data. 


INVERS - LOGICAL 
If INVERS is set to .FALSE. the transform 


ES? 


| | tj —1)(k 1) 
X;+2Y; = A B 21 
j j 2a k ai z)exp( 2m ES? ) 
is calculated, where j = 1, 2,... , ES? and the summation is over k = 1, 2,... ,ES?. 
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If INVERS is set to .TRUE. the transform 


ES? ra 
. i(j — 1)(k — 1) 
Ax + iB; = 2 (X; + #Y; Jeap(—20 5") 
is calculated,where k = 1, 2, ..., ES? and the summation is over j = 1, 2,..., ES?. 


The argument is unchanged on exit. 
FIRST —- LOGICAL 


If FIRST is set to .TRUE. the exponential coefficients for the transform are calculated. 
Consequently FIRST must be set to .TRUE. the first time this routine is called within 
a program,but may be set to .FALSE. for all subsequent calls. 


The argument is unchanged on exit. 


6 Error Indicators 


None. 


7 Auxiliary Routines 


The routine calls the DAP library routines Z-CO6bFFTIDCOEFF, Z-CO6ROWFFT, 
Z.CO6COLFFT, Z-CO6FFTIDBREV. 


8 Accuracy 





Accuracy will be data dependent. You can get some idea of the accuracy by carrying out 
the transform, then carrying out the inverse transform and comparing the results with the 
original data. 


9 Further Comments 


The routine uses a common block with name CCO6FFTLV. Consequently the your program 
must not use a common block with this name. 


10 Keywords 


Fast Fourier Transform 


11 Example 


Not available. 
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5.6 COMP_FFT_2D_8_TO_8 


1 Purpose 


COMP. FFT. 2D. 8. TO_8 calculates the two dimensional fast Fourier transform of an 8-bit 
two’s complement complex image producing an 8-bit two’s complement complex image. The 
input data may be scaled; the output transform is scaled automatically. 


2 Specification 


SUBROUTINE COMP_FFT_2D_8_TO_8 (DATA_IN, DATA_OUT, WORK, 
& NS_SIZE, WE_SIZE, 
& SCALE_FACTOR, INVERS, IFAIL) 


INTEGER*1 DATA_IN (,, NS.SIZE,WE-_SIZE,2) 
& DATA. OUT (,, NS_SIZE,WE-SIZE,2) 


REAL*3 WORK (, , NS_SIZE,WE-SIZE,2) 
INTEGER*4 NS_SIZE, WE.SIZE, SCALE_FACTOR, IFAIL 


LOGICAL INVERS 


3 Description 


The two dimensional transform is calculated by performing independent row and column 
transforms on the image. The forward Fourier transform is defined to be: 


Ni —1No—1 


- —272 ,, : a 
X (ki, ke) = oD > exp {ae (bad + ini} X (ji; Ja) 
ji=0 ja=0 


where: 
Ni = NS_SIZE x ES 
No=WE_SIZE x ES 


If INVERS is .TRUE. then the complex conjugate of the exponential is used. Normalisation 
is not performed for either the forward or the inverse transform. 


4 References 
None. 


5 Arguments 
DATA_IN 


On entry DATA.IN contains the image to be transformed. The dimensions of the image 
must be specified in NS_SIZE and WE_SIZE. The image is assumed to be stored in 
DATA_IN using a crinkled mapping where the third and fourth co-ordinates of DATA_IN 
are the vertical and horizontal co-ordinates of each crinkled tile. The real part of the 
image is stored in the first half of DATA.IN, and the imaginary part in the second half. 


Unchanged on exit, unless DATA_IN has been EQUIVALENCEd to DATA.OUT or 
WORK. 
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DATA. OUT 
On exit DATA_OUT contains the two dimensional Fourier transform of DATA_IN. 


WORK 
Used as work space. 


NS.SIZE 


On entry NS_SIZE contains the height in pixels of the image, divided by ES. NS_SIZE 
must be in the range 1 to 32 and a power of two. Unchanged on exit. 


WE_SIZE 


On entry WE_SIZE contains the width of the image in pixels divided by ES. WE_SIZE 
must be in the range 1 to 32 and a power of two. Unchanged on exit. 


SCALE. FACTOR 
On entry SCALE. FACTOR contains the scale-factor by which the input data has been 
scaled. (positive if multiplied by; negative if divided by). On exit SCALE.FACTOR 
contains the factor by which the result has been scaled. 

INVERS 
If INVERS is set to .TRUE. the inverse discrete Fourier transform is -performed. Un- 
changed on exit. 

IFAIL 


TFAIL equals 0 on exit unless an error occurs. (See 6 below) 





6 Errors 
IFAIL=0 Successful exit 
IFAIL=1 NS.SIZE out of range 
IFAIL=2 WE-SIZE out of range 


7 Auxiliary Routines 


This routine calls the Image Processing Library subroutines COMP_FFT. 2D.REAL.3, 
X05LOG2 and X05PI, and the parallel data transform library is also used. 


8 Accuracy 


Not available. 


9 Further Comments 


The transform is calculated in place; DATA-IN, DATA.OUT and WORK may all occupy 
disjoint areas of DAP store or any two may be EQUIVALENCEd or all three may be EQUIV- 
ALENCEd to start at the same address. For example: 


EQUIVALENCE (DATA_IN,DATA_ OUT) 
EQUIVALENCE (DATA.IN,DATA_OUT,WORK) 


would both work. But 





EQUIVALENCE (DATA. IN, WORK),(DATA. OUT, WORK(,,NS-_ SIZE,3)) 
would not work. 
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10 Keywords 


Fast Fourier Transform, Image Processing 
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5. COMP_FFT_2D_8_TO_16 


1 Purpose | 
COMP_FFT_2D_8.TO_16 calculates the two dimensional fast Fourier transform of an 8-bit 


two’s complement complex image producing a 16-bit two’s complement complex image. The 
input data may be scaled; the output transform is scaled automatically. 
2 Specification 
SUBROUTINE COMP.FFT_2D-.8_TO-_16 (DATA_.IN, DATA_OUT, WORK, 
& : NS_SIZE, WE_SIZE, 
& SCALE_FACTOR, INVERS, IFAIL) 
INTEGER*1 DATA_IN (, , NS_SIZE,WE-SIZE,2) 
INTEGER*2 DATA. OUT (,, NS_SIZE,WE_SIZE,2) 
REAL*3 WORK (,, NS_SIZE,WE- SIZE,2) 
INTEGER*4 NS_SIZE, WE.SIZE, SCALE_FACTOR, IFAIL 


LOGICAL INVERS 


3 Description 





The two dimensional transform is calculated by performing independent row and column 
transforms on the image. The forward Fourier transform is defined to be: 





; Ny=-1No—-1 _ 99s 
X (ki, k2) = yD > exp { NN; (hij + tain) } X (j1, J2) 


ji =0 jqx0 


where: 


Ni =NS_SIZE x ES 

N2=WE_SIZE x ES 
If INVERS is .TRUE. then the complex conjugate of the exponential is used. Normalisation 
is not performed for either the forward or the inverse transform. 


4 References 
None. 


5 Arguments 
DATA_IN 


On entry DATA_IN contains the image to be transformed. The dimensions of the image 
must be specified in NS.SIZE and WE_SIZE. The image is assumed to be-stored in 
DATA.IN using a crinkled mapping where the third and fourth co-ordinates of DATA_IN 
are the vertical and horizontal co-ordinates of each crinkled tile. The real part of the 
image is stored in the first half of DATA_IN, and the imaginary part in the second half. 
Unchanged on exit, unless DATA_IN has been EQUIVALENCEd to DATA_OUT or 
WORK. 
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DATA_ OUT 
On exit DATA. OUT contains the two dimensional Fourier transform of DATA _IN. 


WORK 


Used as work space. 


NS.SIZE 


On entry NS_SIZE contains the height in pixels of the image divided by F.S. NS_SIZE 
must be in the range 1 to 32 and a power of two. Unchanged on exit. 


WE.SIZE 


On entry WE_SIZE contains the width of the image in pixels divided by HS. WE_SIZE 
must be in the range 1 to 32 and a power of two. Unchanged on exit. 


SCALE. FACTOR 


On entry SCALE. FACTOR contains the scale-factor by which the input data has been 
scaled. (positive if multiplied by; negative if divided by). On exit SCALE. FACTOR 
contains the factor by which the result has been scaled. 


INVERS 


If INVERS is set to .TRUE. the inverse discrete Fourier transform is performed. Un- 
changed on exit. . 


IFAIL 


IFAIL equals 0 on exit unless an error occurs. (See 6 below) 


Errors 
IFAIL=0 Successful exit 
IFAIL=1 NS.SIZE out of range 
IFAIL=2 WE-SIZE out of range 


Auxiliary Routines 


This routine calls the Image Processing Library subroutines COMP_FFT_2D_.REAL-_3, 
X05LOG2 and XO5PI, and the parallel data transform library is also used. 


Accuracy 


Not available. 


Further Comments 
The transform is calculated in place; DATA-IN, DATA_OUT and WORK may all occupy 
disjoint areas of DAP store or any two may be EQUIVALENCEd or all three may be EQUIV- 
ALENCEd to start at the same address. For example: 

EQUIVALENCE (DATA-_IN,DATA_OUT) 

EQUIVALENCE (DATA_IN,DATA.OUT,WORK) 
would both work. But 


EQUIVALENCE (DATA_IN,WORK),(DATA. OUT, WORK(,,NS_SIZE,3)) 
would not work. 
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5.8 COMP_FFT_2D_8_TO_24 


1 Purpose 


COMP. FFT_2D_8_.TO-24 calculates the two dimensional fast Fourier transform of an 8-bit 
two’s complement complex image producing a 24-bit two’s complement complex image. The 
input data may be scaled; the output transform is scaled automatically. 


2 Specification 


SUBROUTINE COMP.FFT.2D_8_TO_24 (DATA_IN, DATA-OUT, WORK, 
& NS.SIZE, WE_SIZE, 
& SCALE_FACTOR, INVERS, IFAIL) 


INTEGER*1 DATA_IN (,, NS.SIZE,WE-SIZE,2) 
INTEGER*3 DATA. OUT (,, NS_SIZE,WE_SIZE,2) 

REAL*3 WORK (,,NS.SIZE,WE-SIZE,2) 

INTEGER*4 NS_SIZE, WE.SIZE, SCALE.FACTOR, IFAIL 


LOGICAL INVERS 


3 Description 


The two dimensional transform is calculated by performing independent row and column 
transforms on the image.The forward Fourier transform is defined to be: 


Ny,-1 Ne—1 


a —27t, i fe 
X(k1, ke) = ) ) exp {aah a ini) X(j1, 52) 
ji=0 ja=0 ae 


where: 
Ni = NS_SIZE x ES 
Nog=WE_SIZE x ES 


If INVERS is .TRUE. then the complex conjugate of the exponential is used. Normalisation 
is not performed for either the forward or the inverse transform. 


4 References 


None. 


5 Arguments 
DATA.IN 


On entry DATA_IN contains the image to be transformed. The dimensions of the image 
must be specified in NS.SIZE and WE.SIZE. The image is assumed to be stored in 
DATA_IN using a crinkled mapping where the third and fourth co-ordinates of DATA_IN 
are the vertical and horizontal co-ordinates of each crinkled tile. The real part of the 
image is stored in the first half of DATA_IN, and the imaginary part in the second half. 


Unchanged on exit, unless DATA_IN has been EQUIVALENCEd to DATA_OUT or 
WORK. 
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DATA. OUT 
On exit DATA_OUT contains the two dimensional Fourier transform of DATA_IN. 


WORK 


Used as work space. 


NS_SIZE 


On entry NS_SIZE contains the height in pixels of the image divided by ES. NS_SIZE 
must be in the range 1 to 32 and a power of two. Unchanged on exit. 


WE_SIZE 


On entry WE-_SIZE contains the width of the image in pixels divided by ES. WE. SIZE 
must be in the range 1 to 32 and a power of two. Unchanged on exit. 


SCALE. FACTOR 


On entry SCALE. FACTOR contains the scale-factor by which the input data has been 
scaled. (positive if multiplied by; negative if divided by). On exit SCALE.FACTOR 
contains the factor by which the result has been scaled. 


INVERS 


If INVERS is set to .TRUE. the inverse discrete Fourier transform is performed. Un- 
changed on exit. 


IFAIL 


IF'AIL equals 0 on exit unless an error occurs. (See 6 below) 





6 Errors 
IF AIL=0 Successful exit 
IFAIL=1 NS.SIZE out of range 
IFATIL=2 WE.SIZE out of range 


7 Auxiliary Routines 


This routine calls the Image Processing Library subroutines COMP_FFT. 2D.REAL-.3, 
X05LOG2 and X05PI, and the parallel data transform library is also used. 


8 Accuracy 


Not available. 


9 Further Comments 
The transform is calculated in place; DATA-IN, DATA_OUT and WORK may all occupy 
disjoint areas of DAP store or any two may be EQUIVALENCE¢d or all three may be EQUIV- 
ALENCEd to start at the same address. For example: 
EQUIVALENCE (DATA_IN,DATA_OUT) 
EQUIVALENCE (DATA_IN,DATA_OUT,WORK) 
would both work. But 





EQUIVALENCE (DATA_IN s»WORK), (DATA. OUT, WO RK(,,NS. SIZE,3)) 
would not work. 
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Fast Fourier Transform, Image Processing 
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5.9 COMP_FFT_2D_16_TO_16 


1 Purpose | 
COMP. FFT_2D_16_TO-_16 calculates the two dimensional fast Fourier transform of a 16- 
bit two’s complement complex image producing a 16-bit two’s complement complex image. 
The input data may be scaled; the output transform is scaled automatically. 
2 Specification 
SUBROUTINE COMP_FFT-_2D_16.TO_16(DATA_IN, DATA.OUT, WORK, 
& NS_SIZE, WE_SIZE, 
& SCALE_FACTOR, INVERS, IFAIL) 


INTEGER*2 DATA_IN (, , NS_SIZE,WE.SIZE,2) 
& DATA_OUT(,, NS_SIZE,WE_SIZE,2) 


REAL*3 WORK (,, NS_SIZE,WE-SIZE,2) 
INTEGER*4 NS_SIZE, WE.SIZE, SCALE_FACTOR, IFAIL 


LOGICAL INVERS 


3 Description 





The two dimensional transform is calculated by performing independent row and column 
transforms on the image.The forward Fourier transform is defined to be: 


: Ni-1No-1 oo 
X(k1,k2) = 5° >> exp (hi + iain) X(j1,J2) 


j1=0 jo=0 


where: 


Ni, =NS.SIZE x ES 

No=WE_SIZE x ES 
If INVERS is .TRUE. then the complex conjugate of the exponential is used. Normalisation 
is not performed for either the forward or the inverse transform. 


4 References 


None. 


5 Arguments 
DATA_IN 


On entry DATA_IN contains the image to be transformed. The dimensions of the image 
must be specified in NS.SIZE and WE_SIZE. The image is assumed to be stored in 
DATA_IN using a crinkled mapping where the third and fourth co-ordinates of DATA.IN 
are the vertical and horizontal co-ordinates of each crinkled tile. The real part of the 
image is stored in the first half of DATA_IN, and the imaginary part in the second half. 
Unchanged on exit, unless DATA_IN has been EQUIVALENCEd to DATA.OUT or 
WORK. 
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DATA. OUT 
On exit DATA. OUT contains the two dimensional Fourier transform of DATA _IN. 


WORK 


Used as work space. 


NS_SIZE 


On entry NS_SIZE contains the height in pixels of the image divided by ES. NS_SIZE 
must be in the range 1 to 32 and a power of two. Unchanged on exit. 


WE_SIZE 


On entry WE_SIZE contains the width of the image in pixels divided by BS. WE_SIZE 
must be in the range 1 to 32 and a power of two. Unchanged on exit. 

SCALE. FACTOR 
On entry SCALE. FACTOR contains the scale-factor by which the input data has been 


scaled. (positive if multiplied by; negative if divided by). On exit SCALE.FACTOR 
contains the factor by which the result has been scaled. 


INVERS 


If INVERS is set to .TRUE. the inverse discrete Fourier transform is performed. Un- 
changed on exit. 


IF AIL 


IFAIL equals 0 on exit unless an error occurs. (See 6 below) 


Errors 
IFAIL=0 Successful exit 
IFAIL=1 NS.SIZE out of range 


IFATL=2 WE.SIZE out of range 


Auxiliary Routines 


This routine calls the Image Processing Library subroutines COMP_FFT_.2D.REAL_3, 
X05LOG2 and XO5PI, and the parallel data transform library is also used. 


Accuracy 


Not available. 


Further Comments 
The transform is calculated in place; DATA.IN, DATA.OUT and WORK may all occupy 
disjoint areas of DAP store or any two may be EQUIVALENCEd or all three may be EQUIV- 
ALENCEd to start at the same address. For example: 

EQUIVALENCE (DATA.IN,DATA. OUT) 

EQUIVALENCE (DATA_IN DATA_OUT,WORK) 
would both work. But 


EQUIVALENCE (DATA_IN,WORK),(DATA. OUT,WORK(,,NS_SIZE,3)) 
would not work. 
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Purpose 


COMP. FFT_2D_16. TO. 24 calculates the two dimensional fast Fourier transform of a 16-— 
bit two’s complement complex image producing a 24-bit two’s complement complex image. 
The input data may be scaled; the output transform is scaled automatically. 


Specification 
SUBROUTINE COMP_FFT.2D.16.TO.24(DATA_IN, DATA.OUT, WORK, 
& | NS_SIZE, WE_SIZE, 
& SCALE. FACTOR, INVERS, IFAIL) 


INTEGER*2 DATA_IN (, , NS_SIZE,WE_SIZE,2) 
INTEGER*3 DATA. OUT (,, NS_SIZE,WE_SIZE,2) 

REAL*3 WORK (,,NS_SIZE,WE-SIZE,2) 

INTEGER*4 NS_SIZE, WE_SIZE, SCALE. FACTOR, IFAIL 


LOGICAL INVERS 


Description 


The two dimensional transform is calculated by performing independent row and column 
transforms on the image. The forward Fourier transform is defined to be: 


Ni~1Ng—-1 


is —271,, , — 
X(k1, ke) = ye >> exp {a (had + kxia)} X (ji, Jj2) 
ji=0 jog=0 ae 


where: 
Ni = NS_SIZE x ES 
No=WE_SIZE x ES 


If INVERS is .TRUE. then the complex conjugate of the exponential is used. Normalisation 
is not performed for either the forward or the inverse transform. 


References 


None. 


Arguments 
DATA_IN 


On entry DATA-_IN contains the image to be transformed. The dimensions of the image 
must be specified in NS_SIZE and WE.SIZE. The image is assumed to be stored in 
DATA_IN using a crinkled mapping where the third and fourth co-ordinates of DATA_IN 
are the vertical and horizontal co-ordinates of each crinkled tile. The real part of the 
image is stored in the first half of DATA.IN, and the imaginary part in the second half. 


Unchanged on exit, unless DATA.IN has been EQUIVALENCEd to DATA.OUT or 
WORK. 
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DATA. OUT 
On exit DATA_OUT contains the two dimensional Fourier transform of DATA_IN. 


WORK 


Used as work space. 


NS_SIZE 


On entry NS_SIZE contains the height in pixels of the image divided by ES. NS_SIZE 
must be in the range 1 to 32 and a power of two. Unchanged on exit. 


WE. SIZE 


On entry WE_SIZE contains the width of the image in pixels divided by ES. WE_SIZE 
must be in the range 1 to 32 and a power of two. Unchanged on exit. 


SCALE. FACTOR 


On entry SCALE. FACTOR contains the scale-factor by which the input data has been 
scaled. (positive if multiplied by; negative if divided by). On exit SCALE_FACTOR 
contains the factor by which the result has been scaled. 


INVERS 


If INVERS is set to .TRUE. the inverse discrete Fourier transform is performed. Un- 
changed on exit. 


IFAIL 


IFAIL equals 0 on exit unless an error occurs. (See 6 below) 





6 Errors 
IFAIL=0 Successful exit 
IFAIL=1 NS.SIZE out of range 
IFAIL=2 WE_SIZE out of range 


7 Auxiliary Routines 


This routine calls the Image Processing Library subroutines COMP.FFT_2D. REAL. 3, 
X05LOG2 and XO05PI, and the parallel data transform library is also used. 


8 Accuracy 


Not available. 


9 Further Comments 
The transform is calculated in place; DATA-IN, DATA_OUT and WORK may all occupy 
disjoint areas of DAP store or any two may be EQUIVALENCE or all three may be EQUIV- 
ALENCEd to start at the same address. For example: 
EQUIVALENCE (DATA_IN,DATA.OUT) 
EQUIVALENCE (DATA_IN,DATA. OUT,WORK) 
would both work. But 





EQUIVALENCE (DATA_IN s»WORK),(DATA. OUT,WORK(,,NS_SIZE,3)) 
would not work. 
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10 Keywords 


Fast Fourier Transform, Image Processing 
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5.11 COMP_FFT_2D_24.T0_24 


1 Purpose 
COMP. FFT_2D_24.TO_ 24 calculates the two dimensional fast Fourier transform of a 24- 
bit two’s complement complex image producing a 24-bit two’s complement complex image. 
The input data may be scaled; the output transform is scaled automatically. 
2 Specification 
SUBROUTINE COMP_FFT_2D_24.TO_24(DATA_IN, DATA.OUT, WORK, 
& NS.SIZE, WE. SIZE, 
& : SCALE. FACTOR, INVERS, IFAIL) 


INTEGER*3 DATA.IN (, , NS_SIZE,WE-_SIZE,2) 
& DATA_OUT (, , NS_SIZE,WE-SIZE,2) 


REAL*3 WORK (,,NS_SIZE,WE-SIZE,2) 
INTEGER*4 NS_SIZE, WE_SIZE, SCALE. FACTOR, IFAIL 


LOGICAL INVERS 


3 Description 





The two dimensional transform is calculated by performing independent row and column 
transforms on the image. The forward Fourier transform is defined to be: 


7 Ni ~1 Ng—1 Oni 
X (ki, ke) = ie > cop | (bj + aja) | XG) 
ji=0 ja=0 


where: 


Ni,=NS.SIZE x ES 

No=WE_SIZE x ES 
If INVERS is .TRUE. then the complex conjugate of the exponential is used. Normalisation 
is not performed for either the forward or the inverse transform. 


4 References 


None. 


5 Arguments 
DATA_IN 


On entry DATA_IN contains the image to be transformed. The dimensions of the image 
must be specified in NS.SIZE and WE_SIZE. The image is assumed to be stored in 
DATA-_IN using a crinkled mapping where the third and fourth co-ordinates of DATA_IN 
are the vertical and horizontal co-ordinates of each crinkled tile. The real part of the 
image is stored in the first half of DATA_IN, and the imaginary part in the second half. 
Unchanged on exit, unless DATA_IN has been EQUIVALENCEd to DATA_OUT or 
WORK. 
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DATA. OUT 
On exit DATA. OUT contains the two dimensional Fourier transform of DATA_LIN. 


WORK 


Used as work space. 


NS_SIZE 


On entry NS_SIZE contains the height in pixels of the image divided by ES. NS_SIZE 
must be in the range 1 to 32 and a power of two. Unchanged on exit. 


WE.SIZE 


On entry WE_SIZE contains the width of the image in pixels divided by HS. WE. SIZE 
must be in the range 1 to 32 and a power of two. Unchanged on exit. 


SCALE. FACTOR 


On entry SCALE. FACTOR contains the scale-factor by which the input data has been 
scaled. (positive if multiplied by; negative if divided by). On exit SCALE. FACTOR 
contains the factor by which the result has been scaled. 


INVERS 


If INVERS is set to .TRUE. the inverse discrete Fourier transform is performed. Un- 
changed on exit. 


IF AIL 


IFAIL equals 0 on exit unless an error occurs. (See 6 below) 


Errors 
IFAIL=0 Successful exit 
IFAIL=1 NS_SIZE out of range 
IFAIL=2 WE-SIZE out of range 


Auxiliary Routines 


This routine calls the Image Processing Library subroutines COMP. FFT.2D.REAL_3, 
X05LOG2 and XO5PI, and the parallel data transform library is also used. 


Accuracy 


Not available. 


Further Comments 
The transform is calculated in place; DATA-IN, DATA.OUT and WORK may all occupy 
disjoint areas of DAP store or any two may be EQUIVALENCEd or all three may be EQUIV- 
ALENCEd to start at the same address. For example: 

EQUIVALENCE (DATA_IN,DATA_OUT) 

EQUIVALENCE (DATA_IN,DATA.OUT,WORK) 
would both work. But 


EQUIVALENCE (DATA_IN,WORK),(DATA-OUT,WORK(,,NS_SIZE,3)) 
would not work. 
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10 Keywords 


Fast Fourier Transform, Image Processing 








Image Processing Library man014.03 91 





COMP_FFT.2D_REAL.3 Low Level Image Processing Routines 





5.12 COMP_FFT_2D_REAL_3 


1 Purpose 


COMP_FFT. 2D. REAL-3 calculates the two dimensional fast Fourier transform of a 24-bit 
floating point complex image producing a 24-bit floating-point complex image. The input 
data may be scaled; the output transform is scaled automatically. 


2 Specification 


SUBROUTINE COMP._FFT_2D.REAL-_3 (DATA, NS_SIZE, WE_SIZE, 
& INVERS, IFAIL) 


REAL*3 DATA (,,NS.SIZE,WE_SIZE,2) 
INTEGER*4 NS_SIZE, WE.SIZE, SCALE_FACTOR, IFAIL 


LOGICAL INVERS 


3 Description 


The two dimensional transform is calculated by performing independent row and column 
transforms on the image. The forward Fourier transform is defined to be: 





27% 


N,-1Nq—1 = 
{ae (hd +- ini) X (ji, j2) 


X(ki, ke) = > >» exp 


ji=0 ja=0 


where: 


Ni, =NS.SIZE x ES 
No=WE.SIZE x ES 


If INVERS is .TRUE. then the complex conjugate of the exponential is used. Normalisation 
is not performed for either the forward or the inverse transform. 


4 References 


None. 


§ Arguments 
DATA 


On entry DATA contains the image to be transformed. The dimensions of the image 
must be specified in NS.SIZE and WE_SIZE. The image is assumed to be stored using a 
crinkled mapping where the third and fourth co-ordinates of DATA are the vertical and 
horizontal co-ordinates of each crinkled tile. The real part of the image is stored in the 
first half of DATA, and the imaginary part in the second half. On exit, DATA contains 
the transformed image. 


NS.SIZE 


On entry NS_SIZE contains the height in pixels of the image divided by BS. NS_SIZE 
must be in the range 1 to 32 and a power of two. Unchanged on exit. 
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WE. SIZE 


On entry WE_SIZE contains the width of the image in pixels divided by ES. WE_SIZE 
must be in the range 1 to 32 and a power of two. Unchanged on exit. 


INVERS 


If INVERS is set to TRUE. the inverse discrete Fourier transform is performed. Un- 
changed on exit. 


IF AIL 


IFAIL equals 0 on exit unless an error occurs. (See 6 below) 


6 Errors 
IFAIL=0 Successful exit 
IFAIL=1 NS.SIZE out of range _ 
IFAIL=2 WE.SIZE out of range 


7 Auxiliary Routines 


This routine calls the Image Processing Library subroutines XO5LOG2 and XO5PI, and the 
parallel data transform library is also used. 


8 Accuracy 


Not available. 





9 Further Comments 


None. 


10 Keywords 


Fast Fourier Transform, Image Processing 
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5.13 CONVOLVE.8 


1 Purpose 


CONVOLVE.8 calculates the convolution of an 8-bit per pixel two’s complement image with 
a square mask. The result may be scaled to avoid overflow. 


2 Specification 


SUBROUTINE CONVOLVE8(IMAGE,NS_SIZE, WE_SIZE, CONV, 
& SCALE. FACTOR,MASK_SIZE,MASK,IFAIL) 


INTEGER*1 IMAGE(,,NS.SIZE,WE_SIZE), CONV (,, NS_SIZE,WE_SIZE) 
& MASK (MASK_SIZE,MASK_SIZE) 


INTEGER*4 NS_SIZE, WE.SIZE, MASK_SIZE, SCALE_ FACTOR, IFAIL 


3 Description 


CONVOLVE.8 calculates the convolution of a square user-supplied mask with an image. 
Here, by convolution we mean in the image processing sense, that is a correlation. The 
image is assumed to be extended with the grey-scale value 0 for calculating the convolution. 
Each value returned in CONV is calculated the mask centred over the corresponding pixel. 





4 References 


None. 


5 Arguments 
IMAGE 


On entry IMAGE contains the image to be convolved. Each dimension of the image size 
must be in the range ES to 32HS in multiples of ES . The image is assumed to be stored 
using a crinkled mapping where the third and fourth co-ordinates of IMAGE are the 
vertical and horizontal co-ordinates of each crinkled tile. 


NS.SIZE 
On entry NS_SIZE contains the height in pixels of the image divided by ES. NS_SIZE 
must be in the range 1 to 32. Unchanged on exit. 

WE. SIZE 


On entry WE-_SIZE contains the width of the image divided by ES. WE_SIZE must be 
in the range 1 to 32. Unchanged on exit. 


CONV 
On exit CONV contains the result of convolving IMAGE with MASK divided by 
SCALE.FACTOR. 

SCALE. FACTOR 
On entry SCALE. FACTOR contains the scale-factor which the result is scaled down by 


to avoid overflow. SCALE.FACTOR must be in the range —21° to 2!5 — 1 for CON- 
VOLVE.8 and in the range —2°! to 23! — 1 for CONVOLVE_16. Unchanged on exit. 
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MASK_SIZE 


On entry MASK_SIZE contains the size of the convolving mask. MASK SIZE must be 
one of 3,5,7,...31. Unchanged on exit. 


MASK 


On entry MASK contains the mask that is to be convolved with IMAGE. Unchanged on 
exit. 


IFAIL 


IF AIL equals 0 on exit unless an error occurs. (See 6 below) 


6 Errors 
IFAIL=0 Successful exit 
IFAIL=1 NS_SIZE out of range 
IFAIL=2 WE_SIZE out of range 
IFAIL=3 MASK_SIZE out of range 


IFAIL=4 Overflow in 16-bit to 8-bit conversion 


7 Auxiliary Routines 


The Image Processing Library routines ADD_16, SCAL.MULT_8_16, SCAL-_DIV-. 16, 
SHIFT_COL.NORTH.P, SHIFT. COL.SOUTH.P, SHIFT.COL. WEST-P, 
SHIFT.COL.EAST.P, and CONTRACT-_16.TO-8 are called. 





8 Accuracy 


Not available. 


9 Further Comments 


None. 


10 Keywords 


Convolution, Image Processing 
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1 


Purpose 


DIFF.OF.GAUSS_8 convolves one of a set of nine difference-of-Gaussian masks with an 
8-bit two’s complement image. 


Specification 


SUBROUTINE DIFF.OF_.GAUSS_8 (IMAGE,NS_SIZE, WE-SIZE, DOG, 
& CHANNEL, IFAIL) 


INTEGER*1 IMAGE(,,NS.SIZE,WE_SIZE), DOG(,, NS_SIZE,WE_SIZE) 


INTEGER*4 NS_SIZE, WE.SIZE, CHANNEL, IFAIL 


Description 


DOG is twice the convolution (in the image processing sense) of IMAGE with the selected 
difference-of-Gaussian mask which is specified in CHANNEL. IMAGE is assumed to be 
extended by a field of zeros for the convolution. The difference-of-Gaussian mask is: 





where o takes the following values depending on the value of CHANNEL: 

CHANNEL = 1 2 3 4 5 6 7 8 9 

o = V2 2 2/2 4 4/2 8 8/2 16 16/2 

The mask is taken to be zero for x* + y? > (3c)*; thus a square mask of linear size 60 + 1 


is used. For the large mask sizes IMAGE is scaled down by averaging 2x2, 4x4 and 8x8 
blocks of pixels so that the largest mask size actually used in the computation is 19. 


References 
(1] Marr, D. and Hildreth, E. 
1980 Theory of Edge Detection. 
Proc. R. Soc. Lond. B 207 187-217 
Arguments 
IMAGE 
On entry IMAGE contains the image to be convolved with the difference-of-Gaussian 
mask. Each dimension of the image size must be in the range ES to 32HSin multiples 
of ES. The image is assumed to be stored in IMAGE using a crinkled mapping where 
the third and fourth co-ordinates of IMAGE are the vertical and horizontal co-ordinates 
of each crinkled tile. Unchanged on exit. 
man014.03 AMT 




















Low Level Image Processing Routines | DIFF. OF. GAUSS_ 8 


10 


NS.SIZE 


On entry NS_SIZE contains the height in pixels of the image divided by ES. NS_SIZE 
must be in the range 1 to 32. Unchanged on exit. 


WE.SIZE 


On entry WE-_SIZE contains the width of the image divided by ES. WE_SIZE must be 
in the range 1 to 32. Unchanged on exit. 


DOG 


On exit DOG contains twice the convolution of IMAGE with the difference-of-Gaussian 
mask selected by CHANNEL. DOG stores the convolution using a crinkled mapping (see 
IMAGE above). 


CHANNEL 


On entry CHANNEL contains an integer in the range 1 to 9 which selects the difference- 
of-Gaussian mask which is to be used. If CHANNEL is 4 or 5, then NS_SIZE and 
WE-SIZE must be divisible by 2, if CHANNEL is 6 or 7 NS_SIZE and WE_SIZE must 
be divisible by 4, or if CHANNEL is 8 or 9, NS.SIZE and WE_SIZE must be divisible 
by 8. Unchanged on exit. | 


IFAIL 


IF'AIL has the value 0 on exit unless an error occurs. (See 6 below.) 


Errors 
IFAIL=0 Successful exit 
IFAIL=1 NS_SIZE out of range 
IFAIL=2 WE.SIZE out of range 
IFAIL=6 CHANNEL out of range 
IFAIL=19 NS.SIZE incompatible with CHANNEL 
IFAIL=20 WE.SIZE incompatible with CHANNEL 


Auxiliary Routines 


This routine uses the Image Processing Library routines APPLY. DOG_8, EXPAND_IMAGE.8, 
REDUCE_IMAGE.8 and SHIFT. IMAGE_8. 


Accuracy 


Not available. 


Further Comments 


None. 


Keywords 


Convolution, Difference-of-Gaussians, Image Processing 
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1 Purpose 


F01.G.MM performs a general matrix multiply of two matrices A and B where A is a PxQ 
matrix and B is a QxR matrix with P,Q and R in the range 1 to ES. 


2 Specification 


REAL MATRIX FUNCTION F01-G_MM(A,B,P,Q,R, IFAIL ) 
REAL A(,) , B(,) 


INTEGER P ,Q,R, IFAIL 


3 Description 
The routine is an optimised general matrix multiply using one of the following three proce- 
dures,depending on the relative sizes of P,Q and R (see [1)). 


Procedure 1 


F01GMM=0.0 
DO 10 I=1,Q 
10 FOIGMM=F01GMM+MATC(A(,I))*MATR(B(L,)) 


Procedure 2 
DO 10 I=1,P 
10 FOIGMM(I,)=SUMR(MATC(A(I,))*B) 
Procedure 3 
DO 10 I=1,R 
10 FO1GMM(,I)=SUMC(A*MATR(B(,]))) 


If P/Q > 0.75 and R/Q > 0.75 procedure 1 is used, otherwise if P > R procedure 3 is used 
or if P< R procedure 2 is used; the number 0.75 was determined empirically. 


4 References 


(i] MCKEOWN J.J. 
Multiplication of non-standard matrices on DAP. 
DAP newsletter no. 7. 


5 Arguments 
A ~ REAL MATRIX 


On entry A contains the first of the two matrices to be multiplied together - array 
elements outside the matrix to be multiplied must be set to zero. The contents of A are 
unchanged on exit. 
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B —- REAL MATRIX 


On entry B contains the second of the two matrices to be multiplied together - array 
elements outside the matrix to be multiplied must be set to zero. The contents of B are 
unchanged on exit. 


P ~ INTEGER 
The number of rows in the first matrix. Unchanged on exit. 
Q - INTEGER 


The number of columns in the first matrix and the number of rows in the second matrix. 
Unchanged on exit. 


R - INTEGER 


The number of columns in the second matrix. Unchanged on exit. 


IFAIL - INTEGER 


Unless the routine detects an error (see 6 below) IFAIL contains zero on exit. 


6 Errors 
Errors detected by the routine: 


IFAIL = 1 At least one of P, Q or R is not in the range 1 to ES. 


7 Auxiliary Routines 


None. 





8 Accuracy 


You can expect six significant figures. 


9 Further Comments 


None. 


10 Keywords 
Matrix multiply. 
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Purpose 


FO01. M_INV calculates, in place, the inverse of a given N x N matrix with N in the range 1 


to ES. 


Specification 


SUBROUTINE F01_M_INV( A, N, IFAIL ) 
REAL A(,) 


INTEGER N , IFAIL 


Description 


The matrix is inverted using Gauss-Jordan elimination with full pivoting. 


References 


None. 


Arguments 
A — REAL MATRIX 


On entry A contains the matrix to be inverted,which is assumed to be located in the top 
left of A and array elements outside the input matrix must be set to zero. On exit A 
contains the inverse of that matrix. 


N —- INTEGER 


On entry N must be set to the order of the matrix to be inverted. N is unchanged on 
exit. 


IFAIL —- INTEGER 


Unless the routine detects an error (see 6 below) IFAIL contains zero on exit. 


Errors 


Errors detected by the routine: 


IFAIL=1 N is not in the range 1 to BS. 
IFAIL=2 A pivot element is equal to zero — the matrix is singular 


Auxiliary Routines 


None. 


Accuracy 


You can expect five or six significant figures for well conditioned problems. 
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9 Further Comments 


None. 


10 Keywords 


Matrix inversion, Gauss-Jordan elimination. 
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5.17 FO04_QR_GIVENS_SOLVE 


1 Purpose 


F04. QR. GIVENS_ SOLVE solves the linear system Ax = 6 for x where A is an nxn matrix 
with 2<n<33. The routine may be used to simultaneously solve for up to ES different right 
hand side vectors 6. 


2 Specification 
SUBROUTINE F04.QR.GIVENS_SOLVE(A,X,B,N,NB, IFAIL ) 
INTEGER N , NB, IFAIL 


REAL A(,) , X(,) , B(,) 


3 Description 
The routine factorizes the given nxn matrix A as: 
QA=R 
where Q is an orthogonal matrix and r is upper triangular. 


Givens method of plane rotations is used to annihilate elements of A below the leading 
diagonal until the matrix R remains. This leaves an upper triangular system which is solved 
by back substitution. Row 7 of A is used to annihilate the element in position (i+ 1,j) by 
pre-multiplying A by a matrix of the form: 





Pid = diag (IG-1),UGi41), n-i-1)) 1S I< n-1 


GC 


where Us +1 =— ( Sigs, cps 
3 3 


; , with c? +s? =1 
In the usual serial application, these rotations are applied sequentially, but on the DAP you 
can perform up to > rotations simultaneously [1]. 

4 References 


(1] SAMEH A H and KUCK DJ . 
On stable parallel linear system solvers 
Journal of the Association of Computing Machinery, Vol 25, No 1, pp 81-91. 


5 Arguments 


A —- REAL MATRIX 


On entry, elements Ay; ;) (i= 1, 2,...,.N; 7 = 1,2, ..., N) must be set to the elements 
of the matrix defining the linear system. A is unchanged on exit. 
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X —- REAL MATRIX 


On exit, column 7 of X will contain the solution of the system corresponding to the jth 
column of B. 


B —- REAL MATRIX 


On entry, columns 1 to NB must give the NB right hand side vectors. B is unchanged 
on exit. 


N —- INTEGER 


On entry, N must be set to the order of the matrix A. N is unchanged on exit. 


NB - INTEGER 


On entry, NB must be set to the number of right hand side vectors for which the system 
is to be solved. NB is unchanged on exit. 


IFAIL - INTEGER 


Unless the routine detects an error (see section 6) IFAIL contains zero on exit. 


Errors 
Errors detected by the routine: 
IFAIL = 1 N is not in the range 3 to HSor NB is not in the range 1 to ES 
IFAIL = 2 A zero pivot has been found during the back substitution process, 
| that is, the matrix is singular 
IFAIL = 3 A very small pivot has been found during the back substitution process 
and the matrix is probably singular. Computation proceeds anyway, 
but you should treat the results with caution 
Auxiliary Routines 


This routine calls library routines Z_.F04. BACK. SUBST, 
Z.FO4_SPREAD.LMAT-_ EAST, Z_FO4_.SPREAD.RMAT-_EAST and Z.FO4_ UPDATE. 
Accuracy 
Empirical results indicate that errors may be expected in the gth Gy 7th significant digit. 
The routine will return IFAIL = 3 (see 6 below) if the condition: 

MAXi,; |Rij| 


>5x 105 
MIN; |Rj;| 


is satisfield, where R;; is the upper triangular matrix defined in section 3. 


Further Comments 


You must not use common blocks with the names: 


C.F04.QR1 and C_F04_QR2 


Keywords 


Givens’ Rotation, Linear Equations. 
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Purpose 


FILL.IN_1 performs a region growing process on a logical image starting from input seed 
points. The region growing is constrained by the input image — it can only grow where the 
image is .TRUE.. 


Specification 


SUBROUTINE FILL.IN.1(BLOBS,NS_SIZE, WE_SIZE, SEED, 
& FILLED, PASSES, IFAIL) 


LOGICAL BLOBS (,,NS.SIZE,WE_SIZE), SEED (, ,NS_SIZE,WE_SIZE) 
& FILLED (, , NS. SIZE,WE- SIZE) 


INTEGER*4 NS_SIZE, WE_SIZE, PASSES, IFAIL 


Description 


FILL.IN.1 grows out a .TRUE. region from each . TRUE. point in SEED. The nearest 
neighbours of each . TRUE. point in the growing region are also set to be . TRUE. provided 
the equivalent point in BLOBS is also .TRUE.. This process is iterated until no further 
growth occurs (or PASSES iterations have been made). Thus FILLED is identical to BLOBS 
for all ‘blobs’ with seed points in them marked in SEED (provided no more than PASSES 
iterations are required to do this). 


References 


None. 


Arguments 
BLOBS 


On entry BLOBS contains logical bit-map image. Each dimension of the image size must 
be in the range ES to 32ESin multiples of ES. The image is assumed to be stored in 
IMAGE using a crinkled mapping where the third and fourth co-ordinates of BLOBS are 
the vertical and horizontal co-ordinates of each crinkled tile. Unchanged on exit. 


NS_SIZE 


On entry NS_SIZE contains the height in pixels of the image divided by ES . NS.SIZE 
must be in the range 1 to 32. Unchanged on exit. 


WE-SIZE 


On entry WE-.SIZE contains the width of the image in pixels divided by ES. WE. SIZE 
must be in the range 1 to 32. Unchanged on exit. 


SEED 


On entry SEED contains a logical bit-map of the same size as BLOBS. SEED must be 
set to .TRUE. at each seed point to be grown from. SEED is assumed to be stored using 
a crinkled mapping (see BLOBS above). Unchanged on exit. 
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FILLED 


On exit FILLED is identical to BLOBS for all ‘blobs’ with .TRUE. seeds in them. 
FILLED is stored using a crinkled mapping (see BLOBS above). 


PASSES 


On entry PASSES contains the number of iterations of the growing process allowed. 


Unchanged on exit. 
IFAIL 


IFAIL has the value 0 on exit unless an error occurs. (See 6 below.) 


6 Errors 


IFAIL=0 
IFAIL=1 
IFAIL=2 
IFAIL=8 
IFAIL=10 


7 Auxiliary Routines 


None. 


8 Accuracy 





Not applicable. 


9 Further Comments 


None. 


10 Keywords 


Successful exit 

NS.SIZE out of range 

WE.SIZE out of range 

PASSES not positive 

Area fill incomplete after PASSES iterations 


Area Fill, Image Processing 
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Purpose 


HISTOGRAM. C._8 calculates the histogram for an 8-bit two’s complement crinkle mapped 
image or tiles thereof. 


Specification 
SUBROUTINE HISTOGRAM.C_8 (IMAGE, NS_SIZE, WE_SIZE, HISTS, 
& HIST_SIZE, NO_HISTS, BINS, IFAIL) 
INTEGER*1 IMAGE(,,NS_SIZE,WE_SIZE) 
INTEGER*4 NS_SIZE, WE_SIZE, HISTS (BINS,NO_HISTS), HIST. SIZE, 
& NO.HISTS, BINS, IFAIL 
Description 


HISTS is divided into a number of bins, each of which will contain the number of pixels 
with grey-scale values within its range. The number of bins is specified in BINS, and the 
number of grey-scale values in each range is 256 divided by BINS. For example, if BINS 
is 32, pixels with grey-scale values from —128 to -121 will be counted in BIN), those with 
grey-scale values from —120 to -113 will be counted in BINg, etc. 


The sub-images histogrammed are square tiles with sides of ES x HIST_SIZE. 


References 
None. 


Arguments 
IMAGE 


On entry IMAGE contains the image to be histogrammed. Each dimension of the image 
size must be in the range HE Sto 32ESin multiples of ES. The image is assumed to be 
stored in IMAGE using a crinkled mapping where the third and fourth co-ordinates of 
IMAGE are the vertical and horizontal co-ordinates of each crinkled tile. Unchanged on 
exit. 
NS.SIZE 

On entry NS.SIZE contains the height in pixels of the image divided by ES. NS.SIZE 
must be in the range 1 to 32. If local histograms are required then only powers of two 
are allowed. Unchanged on exit. 


WE.SIZE 


On entry WE-SIZE contains the width of the image in pixels divided by HS. WE.SIZE 
must be in the range 1 to 32. If local histograms are required then only powers of two 
are allowed. Unchanged on exit. 


HISTS 


On exit HISTS contains the histogram results for each sub-image. 
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HIST. SIZE 


On entry HIST_SIZE contains the size of side of the square sub-images to be his- 
togrammed, divided by ES. HIST_SIZE must be a power of two unless a global his- 
togram is required. Unchanged on exit. 


NO.HISTS 


On entry NO.HISTS contains the number of histograms to be calculated. NO. HISTS 
must be equal to (NS_SIZE x WE_SIZE)/HIST.SIZE?, unless NO_HISTS is 1, in 
which case a global histogram is calculated. Unchanged on exit. 


BINS 


On entry BINS is the number of ranges of grey-scale values into which the image is to 
be divided. BINS must be one of 32, 64, 128 or 256 for both HISTOGRAM.C_8 and 
HISTOGRAM.C.16. Unchanged on exit. 


IFAIL 


IFAIL has the value 0 on exit unless an error occurs. (See 6 below) 


Errors 
IFAIL=0 Successful exit 
IFAIL=1 Invalid NS_SIZE 
IFAIL=2 Invalid WE_SIZE 
IFAIL=14 Invalid HIST. SIZE 
IFAIL=15 HIST_SIZE incompatible with NO_HISTS 
IFAIL=16 Invalid BINS 


Auxiliary Routines 


None. 


Accuracy 
Not applicable. 


| Further Comments 


None. 


Keywords 


Histogram, Image Processing 
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5.20 HISTOGRAM_S_8 


1 Purpose 


HISTOGRAM.S.8 calculates the histogram for an 8-bit two’s complement sheet ameped 
image or tiles thereof. 


2 Specification 


SUBROUTINE HISTOGRAM_S_8 (IMAGE, NS_SIZE, WE_SIZE, HISTS, 
& HIST_SIZE, NO_HISTS, BINS, IFAIL) 


INTEGER*1 IMAGE(,, NS.SIZE,WE_SIZE) 


INTEGER*4 NS_SIZE, WE.SIZE, HISTS (BINS,NO_HISTS), HIST. SIZE, 
& NO-.HISTS, BINS, IFAIL 


3 Description 


HISTS is divided into a number of bins, each of which will contain the number of pixels 
with grey-scale values within its range. The number of bins is specified in BINS, and the 
number of grey-scale values in each range is 256 divided by BINS. For example, if BINS 
is 32, pixels with grey-scale values from —128 to —121 will be counted in BIN, those with 
grey-scale values from —120 to -113 will be counted in BINg, etc. 


The sub-images histogrammed are square tiles with sides of ES x HIST_SIZE. 





4 References 


None. 


5 Arguments 
IMAGE 


On entry IMAGE contains the image to be histogrammed. Each dimension of the image 
size must be in the range ESto 32HSin multiples of ES. The image is assumed to 
be stored in IMAGE using a sheet mapping where the third and fourth co-ordinates of 
IMAGE label the sheets in the vertical and horizontal directions. Unchanged on exit. 


NS.SIZE 


On entry NS_SIZE contains the height in pixels of the image divided by ES. NS_SIZE 
must be in the range 1 to 32. Unchanged on exit. 


WE.SIZE 


On entry WE.SIZE contains the width of the image divided by ES. WE_SIZE must be 
in the range 1 to 32. Unchanged on exit. 


HISTS 


On exit HISTS contains the histogram results for each sub-image. 
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HIST SIZE 


On entry HIST.SIZE contains the size of side of the square sub-images to be his- 
togrammed, divided by ES. HIST_SIZE must be a factor of NS_SIZE and WE_SIZE 
unless a global histogram is required. Unchanged on exit. 


NO_HISTS 


On entry NO_HISTS contains the number of histograms to be calculated. NO.HISTS 
must be equal to (NS_SIZE x WE_SIZE)/HIST_SIZE?, unless NO_HISTS is 1, in 
which case a global histogram is calculated. Unchanged on exit. 


BINS 


On entry BINS is the number of ranges of grey-scale values into which the image is to 
be divided. BINS must be one of 32, 64, 128 or 256 for both HISTOGRAM -S.8 and 
HISTOGRAM-_S.16. Unchanged on exit. 


IFAIL 


IFAIL has the value 0 on exit unless an error occurs. (See 6 below.) 





6 Errors 
TFAIL=0 Successful exit 
IFAIL=1 NS_SIZE out of range 
IFAIL=2 WE-SIZE out of range 
IFAIL=14 Invalid HIST_ SIZE 
IFAIL=15 ) NS.SIZE not divisible by HIST. SIZE 
IFAIL=16 WE.SIZE not divisible by HIST_SIZE 
IFAIL=17 HIST.SIZE incompatible with NO. HISTS 
IFAIL=18 Invalid BINS 


7 Auxiliary Routines 
None. 


8 Accuracy 
Not applicable. 


9 Further Comments 


None. 


10 Keywords 
Histogram, Image Processing 
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5.21 KIRSCH_8 


1 Purpose 


KIRSCH-§8 applies a user-selected Kirsch compass gradient mask to an 8-bit per pixel two’s 
complement image. The result may be scaled to avoid overflow. 


2 Specification 


SUBROUTINE KIRSCH. 8(IMAGE,NS_SIZE, WE.SIZE, MASK_NO, 
& KIRSCH, SCALE. FACTOR, IFAIL) 


INTEGER*1 IMAGE(,,NS_SIZE,WE_SIZE), 
& KIRSCH (, ,NS.SIZE,WE-_SIZE) 


INTEGER*4 NS_SIZE, WE_SIZE, MASK_NO, SCALE. FACTOR, IFAIL 


3 Description 


The possible Kirsch compass gradient masks to be convolved (in the image processing sense) 
with IMAGE are: 





MASK.NO=1: 3 3. 8 MASK.NO=2: 3 3. 3 
(North) 3. 0 3 (Northeast) -5 0 3 
-5 -5 —5 ~5 -5 3 

MASK.NO=3: -5 3. 3 MASK.NO=4: -5 -5 3 
(East) -5 OO 3 (Southeast) -5 0 83 
—-5 3 3 3 3 38 

MASK.NO=5: -5 —5 —5 MASK.NO=6: 3 -5 —-5 
(South) 3. 0 3 (Southwest) 3. 0 —5 
3 3 3 3. 3 3 

MASK.NO=7: 3 3 -—-5 MASK_.NO=8: 3 3. 3 
(West) 3 O —5 (Northwest) 3. 0 —5 
3. 3 —5 3 9-5 —5 


4 References 


None. 


5 Arguments 
IMAGE 


On entry IMAGE contains the image to be convolved with the Kirsh compass gradient 
mask. Each dimension of the image size must be in the range ES to 32ESin multiples 
of ES. The image is assumed to be stored in IMAGE using a crinkled mapping where 
the third and fourth co-ordinates of IMAGE are the vertical and horizontal co-ordinates 
of each crinkled tile. Unchanged on exit. 
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NS.SIZE 


On entry NS_SIZE contains the height in pixels of the image divided by ES. NS_SIZE 
must be in the range 1 to 32. Unchanged on exit. 
WE. SIZE : 


On entry WE_SIZE contains the width of the image divided by ES. WE_SIZE must be 
in the range 1 to 32. Unchanged on exit. 


KIRSCH 


On exit KIRSCH contains the convolution of IMAGE with the Kirsch compass gradient 
mask selected by MASK_NO divided by SCALE. FACTOR. 


MASK. NO 


On entry MASK_NO contains an integer in the range 1 to 8 which selects the mask to 
be used. Unchanged on exit. 


SCALE. FACTOR 


On entry SCALE_ FACTOR contains the scale-factor which the result is scaled down by 
to avoid overflow. SCALE_ FACTOR must be in the range —215 to 2!5—1 for KIRSCH_8, 
and in the range —2°* to 23! — 1 for KIRSCH_16. Unchanged on exit. 


IFAIL 


IFAIL equals 0 on exit unless an error occurs. (See 6 below) 


Errors 
IF AIL=0 Successful exit 
IFAIL=1 NS.SIZE out of range 
IFAIL=2 WE-SIZE out of range 
IFAIL=3 MASK_NO out of range 
Auxiliary Routines 


Not available. 


Accuracy 
Not available. 


Further Comments 


None. 


Keywords 


Kirsch, Compass Gradient, Edge Detection, Image Processing 
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Purpose 


LAPLACE. 8 applies a user-selected Laplacian mask to an 8-bit per pixel two’s complement 
image. The result may be scaled to avoid overflow. 


Specification 


SUBROUTINE LAPLACE-. 8 (IMAGE,NS_SIZE, WE.SIZE, MASK_NO, 
& LAP,SCALE_ FACTOR, IFAIL) 


INTEGER*1 IMAGE(,,NS_SIZE,WE_SIZE), LAP (, , NS.SIZE,WE-SIZE) 
INTEGER*4 NS_SIZE, WE.SIZE, MASK_NO, SCALE_ FACTOR, IFAIL 


Description 


The possible Laplacian masks to be convolved (in the image processing sense) with IMAGE 
are: 


MASK.NO=1: 0 —]1 0 MASK_.NO=z2: ~1 ~] —l 
—1 4 —] —1 8 —|l 
i =i 0 7, gl I 
MASK.NO=W73: 1 ~2 1 
2 4 2 
1 —2 1 
References 
None. 
Arguments 
IMAGE 


On entry IMAGE contains the image to be convolved with the Laplacian masks. Each 
dimension of the image size must be in the range HS to 32ESin multiples of ES. The 
image is assumed to be stored using a crinkled mapping where the third and fourth 
co-ordinates of IMAGE are the vertical and horizontal co-ordinates of each crinkled tile. 
Unchanged on exit. 


NS.SIZE 


On entry NS_SIZE contains the height in pixels of the image divided by ES. NS_.SIZE 
must be in the range 1 to 32. Unchanged on exit. 


WE.SIZE 


On entry WE_SIZE contains the width of the image divided by ES. WE.SIZE must be 
in the range 1 to 32. Unchanged on exit. 


LAP 


On exit LAP contains the convolution of IMAGE with the Laplacian mask selected by 
MASK.NO divided by SCALE. FACTOR. 
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MASK.NO 


On entry MASK_NO contains 1, 2, or 3 which selects the Laplacian mask to be used. 
Unchanged on exit. 


SCALE. FACTOR 


On entry SCALE_ FACTOR contains the scale-factor which the result is scaled down by to 
avoid overflow. SCALE. FACTOR must be in the range —215 to 215-1 for LAPLACE-8, 
and in the range —2°! to 23! — 1 for LAPLACE_16. Unchanged on exit. 


IFAIL 


IFAIL equals 0 on exit unless an error occurs. (See 6 below) 


Errors 
TFAIL=0 Successful exit 
IFAIL=1 NS_SIZE out of range 
IFAIL=2 WELSIZE out of range 
IFAIL=3 MASK_NO out of range 
Auxiliary Routines 


Not available. 


Accuracy _ 


Not available. 


Further Comments 


None. 


Keywords 


Laplacian Smoothing, Filtering, Image Processing 
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5.23 LINE._DET_8 


1 Purpose 


LINE. DET_8 applies a user-selected line detection mask to an 8-bit per pixel two’s comple- 
ment image. The result may be scaled to avoid overflow. 


2 Specification 


SUBROUTINE LINE_ DET. 8(IMAGE,NS. SIZE, WE_SIZE, MASK_NO, 
& LINE. DET, SCALE. FACTOR, IFAIL) 


INTEGER*1 IMAGE(,,NS.SIZE,WE-_SIZE), 
& LINE. DET (,, NS.SIZE,WE-SIZE) 


INTEGER*4 NS.SIZE, WE.SIZE, MASK.NO, SCALE. FACTOR, IFAIL 


3 Description 


The possible line detection masks that can be convolved (in the image processing sense) with 





IMAGE are: 
MASK.NO=1: -1 -1 -1 MASK.NO=z=2: -1 —1 2 
2 2 2 —1 2 —1 
—1 -—-1 —l 2-1 —l 
MASK.NO=3: -—1l 2 —l MASK._.NO=4: 2-1 ~—1l 
—1 2 —l —1 2 —1 
—] 2 —l -—1 —1 2 


4. References 


None. 


5 Arguments 
IMAGE 


On entry IMAGE contains the image to be convolved with the line detection mask. Each 
dimension of the image size must be in the range ES to 32ESin multiples of ES. The 
image is assumed to be stored in IMAGE using a crinkled mapping where the third 
and fourth co-ordinates of IMAGE are the vertical and horizontal co-ordinates of each 
crinkled tile. Unchanged on exit. 

NS.SIZE 


On entry NS_SIZE contains the height in pixels of the image divided by ES . NS_SIZE 
must be in the range 1 to 32. Unchanged on exit. 
WE. SIZE 


On entry WE-SIZE contains the width of the image in pixels divided by ES. WE.SIZE 
must be in the range 1 to 32. Unchanged on exit. 
LINE.DET 


On exit LINE. DET contains the convolution of IMAGE with the line detection mask 
selected by MASK.NO divided by SCALE. FACTOR. 
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MASK. NO 


On entry MASK_NO contains an integer in the range 1 to 4 which selects the mask to 
be used. Unchanged on exit. 


SCALE. FACTOR 


On entry SCALE. FACTOR contains the scale-factor which the result is scaled down by to 
avoid overflow. SCALE. FACTOR must be in the range —2!5 to 2!5—1 for LINE_DET_8, 
and in the range —2°! to 23! — 1 for LINE. DET_16. Unchanged on exit. 


IFAIL 


IFAIL equals 0 on exit unless an error occurs. (See 6 below) 


6 Errors 
IFAIL=0 Successful exit 
IFAIL=1 .NS_SIZE out of range 
IFAIL=2 WE.SIZE out of range 
IFAIL=3 MASK_NO out of range 


7 Auxiliary Routines 


Not available. 


8 Accuracy 





Not available. 


9 Further Comments 


None. 


10 Keywords 


Line Detection, Image Processing 
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Purpose 


NORMALIZE_ 8 expands the range of grey-scale values used within an 8-bit two’s comple- 
ment image to the full range. 


Specification 
SUBROUTINE NORMALIZE. 8 ( IMAGE,NS_SIZE, WE_SIZE, NORM, IFAIL) 
INTEGER*1 IMAGE(,, NS_SIZE,WE-SIZE), NORM (, , NS_SIZE,WE- SIZE) 


INTEGER*4 NS.SIZE, WE_SIZE, IFAIL 


Description 


The grey-scale range of IMAGE is expanded by scaling linearly to fill the entire range -128 
to 127. 


References 


- None. 


Arguments 
IMAGE 


On entry IMAGE contains the image to be normalized. Each dimension of the image 
size must be in the range ES to 32ESin multiples of ES. The image is assumed to be 
stored in IMAGE using a crinkled mapping where the third and fourth co-ordinates of 
IMAGE are the vertical and horizontal co-ordinates of each crinkled tile. Unchanged on 
exit. 

NS_SIZE 
On entry NS_SIZE contains the height in pixels of the image divided by BS. NS_SIZE 
must be in the range 1 to 32. Unchanged on exit. 

WE._SIZE 
On entry WE_SIZE contains the width of the image divided by ES. WE_SIZE must be 
in the range 1 to 32. Unchanged on exit. 


NORM 


On exit NORM contains the normalization of IMAGE. NORM stores the convolution 
using a crinkled mapping (see IMAGE above). 


IFAIL 


IFAIL has the value 0 on exit unless an error occurs. (See 6 below) 


Errors 
IFAIL=0 Successful exit 
IFAIL=1 NS_SIZE out of range 
IFAIL=2 WE.SIZE out of range 
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7 Auxiliary Routines 


None. 


8 Accuracy 
Not available. 


9 Further Comments 


None. 


10 Keywords 


Normalization, Image Processing 
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5.25 PERC_THRESH-_8 


1 Purpose 
PERC. THRESH-8 performs thresholding on an 8-bit two’s complement image. 


2 Specification 


SUBROUTINE PERC. THRESH-8(IMAGE, NS_SIZE, WE_SIZE, P.THRESH, 
& HIGH_TIDE, PERCENTILE, FLIP, IFAIL) 


INTEGER*1 IMAGE(, , NS_SIZE,WE-SIZE) 


INTEGER*4 NS_SIZE, WE_SIZE, HIGH. TIDE, PERCENTILE, IFAIL 


LOGICAL P. THRESH (,, NS_SIZE,WE_SIZE), FLIP 


3 Description 


P_THRESH is a logical bit-map of the image which has .TRUE. values where the IMAGE 
grey-scale values exceed HIGH. TIDE and .FALSE. otherwise. (This is reversed if FLIP is 
.FALSE. instead of .TRUE..) HIGH.TIDE is determined by the requirement that PER- 
CENTILE per cent of the values in P. THRESH must be .TRUE.. | 


4 References | 


None. 


5 Arguments 
IMAGE 


On entry IMAGE contains the image to be thresholded. Each dimension of the image 
size must be in the range ES to 32ESin multiples of HS. The image is assumed to be 
stored in IMAGE using a crinkled mapping where the third and fourth co-ordinates of 
IMAGE are the vertical and horizontal co-ordinates of each crinkled tile. Unchanged on 
exit. 
NS.SIZE 
On entry NS_SIZE contains the height in pixels of the image divided by ES. NS.SIZE 
must be in the range 1 to 32. Unchanged on exit. 
WE- SIZE 
On entry WE_SIZE contains the width of the image in pixels divided by HS. WE.SIZE 
must be in the range 1 to 32. Unchanged on exit. 


P. THRESH 


On exit P. THRESH contains a logical bit-map of from the thresholding of IMAGE. 
P. THRESH stores the bit-map using a crinkled mapping (see IMAGE above). 


HIGH. TIDE 
On exit HIGH. TIDE is the threshold with which IMAGE is compared. 
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PERCENTILE 


HIGH.TIDE is adjusted so PERCENTILE percent of the pixel grey-scale values exceed 
HIGH. TIDE (for FLIP .TRUE.). 


FLIP 


If FLIP is .TRUE., P. THRESH is .TRUE. when IMAGE is greater than HIGH. TIDE. 
If FLIP is FALSE., P. THRESH is .FALSE. when IMAGE is greater than HIGH_TIDE. 


IFAIL 


IF'AIL has the value 0 on exit unless an error occurs. (See 6 below.) 


6 Errors 
IFAIL=0 Successful exit 
IFAIL=1 NS_SIZE out of range 
IFAIL=2 WE.SIZE out of range 
IFAIL=17 PERCENTILE out of range 0 to 100 


7 Auxiliary Routines 
ABS. THRESH_8 from the Image Processing Library is called. 


8 Accuracy 
Not applicable. 





9 Further Comments 


None. 


10 Keywords 


Thresholding, Image Processing 
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5.26 PREWITT_8 


1 Purpose 


PREWITTS applies a user-selected Prewitt compass gradient mask to an 8-bit per pixel 
two’s complement image. The result may be scaled to avoid overflow. 


2 Specification 


SUBROUTINE PREWITT-.8 (IMAGE,NS_SIZE, WE_SIZE, MASK_NO, 
& PREW, SCALE. FACTOR, IFAIL) 


INTEGER*1 IMAGE(,,NS_SIZE,WE_SIZE), PREW (,, NS_.SIZE,WE-_SIZE) 


INTEGER*4 NS_SIZE, WE_SIZE, MASK.NO, SCALE. FACTOR, IFAIL 


3 Description 


The possible Prewitt compass gradient masks that can be convolved (in the image processing 
sense) to IMAGE are: 


MASK.NO=1: 1 1 1 MASK.NO=2: 1 1 1 





(North) 1-2 1 (Northeast) ~1 -2 1 
—1 +1 —1. —-1 -l 1 

MASK.NO=3: -1 1 1 MASK.NO=4: -1 -1 1 
(East) -1 -2 1 (Southeast) ~1 -2 1 
-l 21 1 1 1 1 

MASK.NO=5: ~—1 -1 -1 MASK.NO=6: 1 -1 -1 
(South) 1-2 1 (Southwest) 1 —-2 -1 
1 1 1 1 1 1 

MASK.NO=7T: 1 1-1 MASK.NO=8: 1 1 1 
(West) 1 -—2 —-1 (Northwest) 1 -2 -1 
1 1 —1 1 —1 —1 


4 References 
None. 


5 Arguments 


IMAGE 


On entry IMAGE contains the image to be convolved with the Prewitt compass gradient 
mask. Each dimension of the image size must be in the range HSto 32ESin multiples 
of ES. The image is assumed to be stored in IMAGE using a crinkled mapping where 
the third and fourth co-ordinates of IMAGE are the vertical and horizontal co-ordinates 
of each crinkled tile. Unchanged on exit. 
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NS.SIZE 


On entry NS_SIZE contains the height in pixels of the image divided by B.S. NS_SIZE 
must be in the range 1 to 32. Unchanged on exit. 
WE. SIZE 


On entry WE_SIZE contains the width of the image divided by ES. WE_SIZE must be 
in the range 1 to 32. Unchanged on exit. 


PREW 


On exit PREW contains the convolution of IMAGE with the Prewitt compass gradient 
mask selected by MASK_NO divided by SCALE. FACTOR. 


MASK. NO 


On entry MASK.NO contains an integer in the range 1 to 8 which selects the mask to 
be used. Unchanged on exit. 


SCALE. FACTOR 


On entry SCALE_ FACTOR contains the scale-factor which the result is scaled down by to 
avoid overflow. SCALE. FACTOR must be in the range —2!5 to 2!5—1 for PREWITT-_8, 
and in the range —2°* to 23! — 1 for PREWITT_16. Unchanged on exit. 


IFAIL 


IFAIL equals 0 on exit unless an error occurs. (See 6 below) 


Errors 
IFATL=0 ” Successful exit 
IFAIL=1 NS.SIZE out of range 
TFAIL=2 WE.SIZE out of range 
IFAIL=3 MASK_NO out of range 


Auxiliary Routines 


Not available. 


Accuracy 


Not available. 


Further Comments 


None. 


Keywords 


Prewitt, Compass Gradient, Image Processing 
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Purpose 


PSEUDO. MEDIAN.S8 calculates the pseudo-median at each pixel in an 8-bit two’s comple- 
ment image for a rectangular neighbourhood centred on the image. 


Specification 
SUBROUTINE PSEUDO. MEDIAN. 8 (IMAGE,NS_SIZE, WE.SIZE, PS.MED, 
& NS_MED, WE.MED, IFAIL) 
INTEGER*1 IMAGE(,,NS.SIZE,WE_SIZE), 
& PS. MED (,, NS_SIZE,WE_SIZE) 
INTEGER*4 NS_SIZE, WE.SIZE, NS.MED, WE. MED, IFAIL 
Description 


The pseudo-median is calculated by first calculating the median of the rows or columns of 
the rectangular neighbourhood about each point defined by NS.MED and WE_MED and 
then calculating the median of the results. The longer dimension is calculated first. 


References 


None. 


Arguments 
IMAGE 


On entry IMAGE contains the image for which the local pseudo-medians are to be cal- 
culated. Each dimension of the image size must be in the range E'S to 32E/Sin multiples 
of ES. The image is assumed to be stored in IMAGE using a crinkled mapping where 
the third and fourth co-ordinates of IMAGE are the vertical and horizontal co-ordinates 
of each crinkled tile. Unchanged on exit. 


NS_SIZE 


On entry NS_SIZE contains the height in pixels of the image divided by ES. NS_SIZE 
must be in the range 1 to 32. Unchanged on exit. 


WE-_SIZE 


On entry WE-SIZE contains the width of the image divided by ES. WE.SIZE must be 
in the range 1 to 32. Unchanged on exit. 


PS. MED 


On exit PS. MED contains the pseudo-median for a rectangular region about each pixel. 
PS.MED is stored using a crinkled mapping (see IMAGE above). 


NS.MED 


On entry NS.MED contains the height of the pseudo-median rectangle which must be 
one of 1, 3, 7, 15, 31. Unchanged on exit. 
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WE. MED 


On entry WE_MED contains the width of the pe acewgeaed rectangle which must be 
one of 1, 3, 7, 15, 31. Unchanged on exit. 


IFAIL 


IFAIL has the value 0 on exit unless an error occurs. (See 6 below) 


6 Errors 
IFAIL=0 ~ Successful exit 
IFAIL=1 NS_SIZE out of range 
IFAIL=2 WE.SIZE out of range 
IFAIL=11 Invalid NS.MED 
IFAIL=12 Invalid WE._.MED 


7 Auxiliary Routines 
The Image Processing Library routines PURE. MEDIAN-.8, SHIFT_IMAGE-8 are called. 


8 Accuracy 
Not applicable. 


9 Further Comments 


None. 





10 Keywords 


Pseudo-median, Median, Image Processing 
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5.28 PURE. MEDIAN_8 


1 Purpose 


PURE. MEDIAN.8 calculates the true median at each pixel in an 8-bit two’s complement 
image for a rectangular neighbourhood centred on the image. 


2 Specification 


SUBROUTINE PURE-MEDIAN-8 (IMAGE,NS_SIZE, WE_SIZE, MED, 
& NS.MED, WE. MED, IFAIL) 


INTEGER*1 IMAGE(,,NS.SIZE,WE_SIZE), 
& MED(,,NS_SIZE,WE-_SIZE) 


INTEGER*4 NS_SIZE, WE.SIZE, NS.MED, WE_MED, IFAIL 


3 Description 


MED is the median grey-scale value of the pixels in the rectangular neighbourhood about 
each point defined by NS.MED and WE. MED. 


4 References 


None. 





5 Arguments 
IMAGE 


On entry IMAGE contains the image for which the local medians are to be calculated. 
Each dimension of the image size must be in the range ES to 32ES in multiples of HS. 
The image is assumed to be stored in IMAGE using a crinkled mapping where the third 
and fourth co-ordinates of IMAGE are the vertical and horizontal co-ordinates of each 
crinkled tile. Unchanged on exit. 

NS. SIZE 
On entry NS_SIZE contains the height in pixels of the image divided by ES. NS_SIZE 
must be in the range 1 to 32. Unchanged on exit. 

WE. SIZE 


On entry WE-SIZE contains the width of the image divided by ES. WE_SIZE must be 
in the range 1 to 32. Unchanged on exit. 


MED 


On exit MED contains the median of a rectangular region about each pixel. MED is 
stored using a crinkled mapping (see IMAGE above). 


NS.MED 


On entry NS.MED contains the height of the median rectangle which must be one of 1, 
3, 7, 15, 31. Unchanged on exit. 
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WE.MED 


On entry WE_MED contains the width of the median rectangle which must be one of 1, 
3, 7, 15, 31. Unchanged on exit. 


IFAIL 


IFAIL has the value 0 on exit unless an error occurs. (See 6 below.) 


6 Errors 
IFAIL=0 Successful exit 
IFAIL=1 NS_SIZE out of range 
IFAIL=2 WE_SIZE out of range 
IFAIL=11 Invalid NS_MED 
IFAIL=12 Invalid WE. MED 


7 Auxiliary Routines 
The Image Processing Library routine SHIFT_IMAGE_8 is called. 


8 Accuracy 
Not applicable. 


9 Further Comments | | | 
15 x 31, 31 x 15 and 31 x 31 medians cannot be performed by PURE.MEDIAN.16. 


10 Keywords 


Median, Image Processing 
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5.29 ROBERTS_8 


1 Purpose 


ROBERTS.8 applies a user-selected Roberts edge detection mask to an 8-bit per pixel two’s 
complement image. The result may be scaled to avoid overflow. 


2 Specification 


SUBROUTINE ROBERTS.8 (IMAGE,NS_SIZE, WE_SIZE, MASK.NO, 
& ROB, SCALE_FACTOR, IFAIL) 


INTEGER*1 IMAGE(,, NS.SIZE,WE_SIZE), ROB(,, NS_SIZE,WE- SIZE) 


INTEGER*4 NS_SIZE, WE_SIZE, MASK_NO, SCALE. FACTOR, IFAIL 


3 Description 


The possible Roberts masks that can be convolved (in the image processing sense) with 
IMAGE are: 


MASK.NO=1: 1 0 MASK.NO=2: 0 1 
0 -1 -1 0 


The masks are extended by zeros to the left and up and then applied as centred 3x3 masks. 


4 References 


None. 


5 Arguments 
IMAGE 


On entry IMAGE contains the image to be convolved with the Roberts edge detection 
mask. Each dimension of the image size must be in the range EHS to 32ESin multiples 
of ES. The image is assumed to be stored in IMAGE using a crinkled mapping where 
the third and fourth co-ordinates of IMAGE are the vertical and horizontal co-ordinates 
of each crinkled tile. Unchanged on exit. 


NS.SIZE 


On entry NS_SIZE contains the height in pixels of the image divided by ES . NS.SIZE 
must be in the range 1 to 32, Unchanged on exit. 
WE-SIZE 


On entry WE-_SIZE contains the width of the image divided by ES. WE_SIZE must be 
in the range 1 to 32. Unchanged on exit. 


ROB 

On exit ROB contains the convolution of IMAGE with the Roberts 

edge detection mask selected by MASK_NO divided by SCALE. FACTOR. 
MASK.NO 


On entry MASK.NO contains an integer in the range 1 to 2 which selects the mask to 
be used. Unchanged on exit. 
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SCALE. FACTOR 


On entry SCALE. FACTOR contains the scale-factor which the result is scaled down by to 
avoid overflow. SCALE. FACTOR must be in the range —2!° to 215—1 for ROBERTS-_8, 
and in the range —2°' to 28! — 1 for ROBERTS_16. Unchanged on exit. 


IFAIL 


IFAIL equals 0 on exit unless an error occurs. (See 6 below) 


Errors 


IFATL=0 
IFAIL=1 
IFAIL=2 
IFAIL=3 


Auxiliary Routines 


Not available. 


Accuracy 


Not available. 


Further Comments 


None. 


Keywords 


Successful exit 
NS_SIZE out of range 
WE.SIZE out of range 
MASK_NO out of range 


Roberts, Edge Detection, Image processing 


Image Processing Library 
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Purpose 


SOBEL-.8 applies a user-selected Sobel edge detection mask to an 8-bit per pixel two’s 
complement image. If MASK.NO= 1 the horizontal edge detection mask is used and if 
MASK_NO= 2 the vertical edge detection mask is used. The result may be scaled to avoid 
overflow. 


Specification 


SUBROUTINE SOBEL. 8(IMAGE, NS_SIZE, WE. SIZE, MASK_NO, SOB, 
& SCALE_ FACTOR, IFAIL) | 


INTEGER*1 IMAGE(,,NS_SIZE,WE_SIZE), SOB(, , NS_SIZE,WE_SIZE) 


INTEGER NS_SIZE, WE_SIZE, MASK.NO, SCALE. FACTOR, IFAIL 


Description 


The vertical and horizontal masks convolved (in the image processing sense) with image are: 


Horizontal: 1 0 —1 Vertical : 1 2 J 
2 0 =2 | 0 0 0 
10 -l . —~l] —2 —-l 
References 
None. 
Arguments 
IMAGE 


On entry IMAGE contains the image to be convolved with the edge detection masks. 
The image is assumed to be stored using a crinkled mapping where the third and fourth 
co-ordinates of IMAGE are the vertical and horizontal co-ordinates of each crinkled tile. 
Unchanged on exit. 


NS.SIZE 
On entry NS_SIZE contains the height of the image in pixels divided by ES. NS_SIZE 
must be in the range 1 to 32. Unchanged on exit. 

WE. SIZE 
On entry WE-.SIZE contains the width of the image in pixels divided by ES. WE_SIZE 
must be in the range 1 to 32. Unchanged on exit. 

MASK.NO 
On entry MASK_NO, contains 1 if the horizontal edge detection mask is to be used and 
MASK.NO, contains 2 if the vertical edge detection mask is to be used. Unchanged on 
exit. 

SOB 


On exit SOB contains the the convolution of IMAGE with either the vertical or horizontal 
Sobel edge detection mask divided by SCALE_.FACTOR. 
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SCALE. FACTOR 


On entry SCALE_ FACTOR contains the scale-factor which the result is scaled down by 
to avoid overflow. SCALE_FACTOR must be in the range —2)5 to 2!5—1 for SOBEL.8, 
and in the range —2°! to 291 — 1 for SOBEL_16. Unchanged on exit. 


IFAIL 


IFAIL equals 0 on exit unless an error occurs. (See 6 below) 


Errors 
IFAIL=0 Successful exit 
IFAIL=1 NS.SIZE out of range 
IFAIL=2 WE-SIZE out of range 
IFAIL=3 MASK._NO not in range 
Auxiliary Routines 


Not available. 


Accuracy 


Not available. 


Further Comments 


None. 


Keywords 


Sobel, Edge Detection, Image Processing 
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5.31 ZERO_X_8 


1 Purpose 


ZERO. X_8 finds the zero crossings in an 8-bit two’s complement image. 
2 Specification 
SUBROUTINE ZERO_X_8 (IMAGE,NS_SIZE, WE_SIZE, ZERO, IFAIL) 
INTEGER*1 IMAGE(,,NS.SIZE,WE-SIZE) 
INTEGER*4 NS_SIZE, WE_SIZE, IFAIL 
LOGICAL ZERO (,,NS_SIZE,WE-_SIZE) 


3 Description 


ZERO is a logical bit-map of the image which has . TRUE. values where the IMAGE changes 
between positive and negative values. - 


4 References 


None. 





5 Arguments 
IMAGE 


On entry IMAGE contains the image to be searched for zero crossings. Each dimension 
of the image size must be in the range ESto 32ESin multiples of ES. The image is 
assumed to be stored in IMAGE using a crinkled mapping where the third and fourth 
co-ordinates of IMAGE are the vertical and horizontal co-ordinates of each crinkled tile. 
Unchanged on exit. 


NS.SIZE 


On entry NS_SIZE contains the height of the image in pixels divided by ES. NS_SIZE 
must be in the range 1 to 32. Unchanged on exit. 


WE. SIZE 


On entry WE-SIZE contains the width of the image divided by ES. WE_SIZE must be 
in the range 1 to 32. Unchanged on exit. 


ZERO 


On exit ZERO contains a logical bit-map of zero crossings in IMAGE. ZERO stores the 
bit-map using a crinkled mapping (see IMAGE above). 


IFAIL 


IFAIL has the value 0 on exit unless an error occurs. (See 6 below.) 
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6 Errors 
IFAIL=0 Successful exit 
IFAIL=1 NS.SIZE out of range 
TFATL=2 WE_SIZE out of range 
IFAIL=7 No zero crossings found 


7 Auxiliary Routines 


None. 


8 Accuracy 
Not applicable. 


9 Further Comments 


None. 


10 Keywords 


Zero Crossing, Image Processing 
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Chapter 6 


Image Analysis Routines 


These image analysis routines allow the interesting features on an image to found and identified 
automatically. SEGMENT-16 and FEATURES. 16 are available for analysing 16-bit images. 
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Purpose 


LABEL. 16 labels the distinct ‘blobs’ in a binary image, returning a 16-bit two’s complement 
image with all the pixels in a blob having the same value, but each blob is given a different 
value. ; 


Specification 


SUBROUTINE LABEL- 16 (BIN-IMAGE, NS_SIZE, WE-SIZE, LAB_IMAGE, 
& IFAIL) 


LOGICAL BIN. IMAGE(,, NS.SIZE,WE-_SIZE) 
INTEGER*2 LAB_IMAGE(, , NS_SIZE,WE-SIZE) 


INTEGER*4 NS_SIZE,, WE_SIZE, IFAIL 


Description 


A unique value is assigned to each . TRUE. pixel in BIN.IMAGE. These values are copy- 
propagated to . TRUE. nearest neighbours. When two or more values are assigned to a pixel 
it is the highest value which is assigned to that pixel. 


References 


None. 


Arguments 


BIN. IMAGE 


On entry IMAGE contains the image on which the blob labelling is to be performed. 
Each dimension of the image size must be in the range 2S to 32E;S in multiples of BS. 
The image is assumed to be stored in IMAGE using a sheet mapping where the third 
and fourth co-ordinates of IMAGE are the vertical and horizontal co-ordinates of each 
sheet. Unchanged on exit. 


NS.SIZE 


On entry NS_SIZE contains the height of the image in pixels divided by ES. NS_SIZE 
must be in the range 1 to ES. Unchanged on exit. 


WE-SIZE 


On entry WE-.SIZE contains the width of the image divided by ES. WE_SIZE must be 
in the range 1 to ES. Unchanged on exit. 


LAB_IMAGE 


On exit LAB_IMAGE contains a sheet-mapped (see IMAGE above) 16-bit image la- 
belling the blobs in BIN.IMAGE. 


IFAIL 


IFAIL has the value 0 on exit unless an error occurs. (See 6 below) 
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6 Errors 
IFAIL=0 Successful exit 
IFAIL=1 NS.SIZE out of range 
IFAIL=2 WE.SIZE out of range. 


7 Auxiliary Routines 
The Image Processing Library routine LABEL_16_SUB is called. 


8 Accuracy 
Not applicable. 


9 Further Comments 


None. 


10 Keywords 


Image Feature Generation, Labelling, Image Processing 
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7 Auxiliary Routines 


None. 


8 Accuracy 
Not available. 


9 Further Comments 


None. 


10 Keywords 


Image Feature Generation, Image Processing 


Table 1 


The following formulae are extracted from Digital Image Processing by Gonzalez and Wintz. 


1. Invariant moment number 2 
eo 2 2 
im2 = (N20 — Moz)” + 401; 
2. Invariant moment number 5 


ims = (N30 = 3712)(N30 + 712) { (nso + m2)* — 3(no1 + nos)" } 
+(3n21 — No3s)(n21 + N03) {3(30 + m2)” — (N21 + N03)" } 


3. Invariant moment number 6 


ime = (n20—mo2){(N30 + M12)" — (21 + N03)"} 
+4711(ns0 + 12)(N21 + Nos) 


4. Invariant moment number 7 


imz = (3n21 —130)(n30 + m2){(n30 + m2)? — 3(na1 + N03)? } 
+(312 — 30)(n21 + 703){3(ns0 + m2)” — (N21 + N03)” } 
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Where: 


a Ung 
PI ™ g(p+q+2)/2 


and a = upp so that the np, are invariant to scale, translation and intensity. 
And where: | 


Upg = ¥_[(z- 2)?(y — 9)" f(z, y)|dady 
f(x,y) is the grey-scale value at position x, y. 


Z,¥y is the position of the centroid of the blob, as defined by: 


_ Yef(z,y)dedy 








r= 
df (x, y)drdy 
es DL Ys (x, y)dzdy 
2. f (2, y)dady 
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